NPS67-82-009 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


j'L*' 


DUAL  PROBE  DIGITAL  SAMPLING  (DPDS) 
TECHNIQUE  FOR  ROTOR  FLOWS —REVIEW 
OF  THE  DEVELOPMENT  AND  AN  ASSESSMENT 
BASED  ON  EARLY  RESULTS 


H.-J.  Heinemann 

♦ 

April  1982 

Final  Report;  October  1980-March  1982 


Approved  for  public  release;  distribution 

unlimited 

/ 

Prepared  for: 

Naval  Air  Systems  Command 

Washington,  DC  20261 

83 

01 

26 

084 

NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  J.  J.  Ekelund  David  A.  Schrady 

Superintendent  Acting  Provost 


The  work  reported  herein  was  carried  out  during  the  period 
October  1980  to  March  1982.  The  author  was  supported  for  twelve 
months  by  the  DFVLR  (German  Aerospace  Research  Establishment) 
and  by  the  project  under  Contract  No.  N62271-80-M-2214.  The 
work  was  costpleted  while  the  author  was  Naval  Air  Systems  Cosi- 
mand  Visiting  Research  Professor  in  Aeronautics  at  the  Naval 
Postgraduate  School  from  October  1981  to  March  1982;  support 
for  this  period  was  under  Work  Request  N00019-82-WR-21036  as 
part  of  Energy  Utilization  and  Airbreathing  Propulsion,  AIRTASK 
(A3103100/186A) /1R024-03-001,  under  the  cognizance  of  Dr.  G. 
Heiche,  AIR  310. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 

This  report  was  prepared  by: 


U  tSRCUUu.  Vl/U2JlMGoc*u~ 
Ig-JdXCHltfHEl 
NAVAIR  Visiting  Research 
Professor  in  Aeronautics 


Reviewed  by: 


UNCLASSIFIED 


sectmrrv  classification  of  this  p»  ,e  »?»an  d««  EmtM 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

U  report  numbe*  |2.  oovt  accession  no. 

NPS67-82-009 

3.  *CCI»1ENT*S  CATALOG  NUMIER 

4.  TITLE  fand  SufcMIlal 

Dual  Probe  Digital  Sampling  (DPDS)  Tech¬ 
nique  for  Rotor  Flows — Review  of  the 
Development  and  an  Assessment  Based  on 
Early  Results 

S  TYPE  OF  REPORT  b  PERI 00  COVERED 

iFinal  Report 

October  1980  -March  1982 

S.  PERFORMING  oro.  report  number 

7.  AuTmON<4> 

H.-J.  Heinemann 

t.  CONTRACT  OR  GRAN*  NUMBER'*) 

»  performing  Organization  NAME  and  ADDRESS 

Naval  Postgraduate  School 

Monterey,  California  93940 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  «  WORK  UNIT  NUMBERS 

NO 00 19 8 2WR- 210 36 

<1.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Naval  Air  Systems  Command 

Washington,  D.C.  20361 

11.  report  OATS 

April  1982 

IS.  NUMBER  of  PAOES 

153 

14.  monitoring  AGENCY  NAME  «  AOORESVII  dllterent  Item  Controlling  Office) 

Naval  Postgraduate  School 

Monterey,  California  93940 

IS.  SECURITY  CLASS,  (ol  eh«*  re pen) 

Unclassified 

<S*.  DECLASSIFICATION/ OOWNGNAOING  1 

SCHEDULE  I 

DISTRIBUTION  STATEMENT  fol  Iflla  Report) 

Approved  for  public  release;  distribution  unlimited 

...  j 

17  DISTRIBUTION  STATEMENT  fat  the  abatract  entered  In  Block  30,  It  dllterent  Item  Report)  | 

! 

*9  <EY  WORDS  (Continue  on  r«v«N*  old#  1/  n«c«««crr  «td  Itfwiliff  fry  bloc*  nustbnr) 

,  Compressor  Flow  Fields 
.  Compressor  Measurements 
|  Probe  Techniques 


T  ASS'*  AC  *  'Continue  an  rarer  ee  altfa  It  neceeeary  and  Identity  by  black  number) 

i  ^>The  state  of  the  art  in  the  development  of  the  DPDS  technique 
at  the  Turbopropulsion  Laboratory  at  the  Naval  Postgraduate 
-  School  is  reported.  In  principle,  this  method  allows  the  de¬ 
termination  of  three  dimensional  flow  fields  behind  rotors.  A 
review  of  the  instrumentation  and  hardware  is  given  and  two 
1  different  data  reduction  procedures  developed  earlier  are  de- 
I  scribed.  The  results  of  reducing  one  set  of  rotor  exit  data 
i _ I _ _ 


DO  1473 


EDITiON  0*  t  NOV  SS  IS  OBSOLETE 
'  '  '■  iO  2-n  14-  440  I  ]_ 


& 


SECURITY  CLASSIFICATION  OF  this  PAGE  '»'>>•*»  Date  fnteced) 


_ UNCLASSIFIED _ 

...unity  CLASSIFICATION  OF  Thi5  PAOCfT1*"  D*'m  Bttffd) 


using  the  two  methods  are  compared.  In  order  to  verify  the 
experimental  procedures  at  low  cost  a  bench-top  periodic  flow 
generator  was  used  first  to  provide  signals  similar  in  type 
and  frequency  to  those  from  the  Transonic  Compressor  in  which 
the  rotor  exit  flow  was  measured.  The  data  measured  in  the 
compressor  were  analyzed,  and  problems  encountered  in  the 
routine  application  of  the  reduction  methods  were  identified. 

The  possible  application  of  similar  techniques  in  the  Ro¬ 
tating  Annular  Cascade  (RGG)  at  the  DFVLR  Gottingen  is  as¬ 
sessed  and  recommendations  are  made  for  the  coordination  of 
future  experimental  programs. 


UNCLASSIFIED 


SBCUNITV  CLASSIFICATION  OF  this  FAOCrKTiAn  Dm* 


ABSTRACT 


The  state  of  the  art  in  the  development  of  the  DPOS  tech¬ 
nique  at  the  Turbopropulsion  Laboratory  at  the  Naval  Post¬ 
graduate  School  is  reported.  In  principle,  this  method  allows 
the  determination  of  three  dimensional  flow  fields  behind 
rotors.  A  review  of  the  instrumentation  and  hardware  is  given 
and  two  different  data  reduction  procedures  developed  earlier 
are  described.  The  results  of  reducing  one  set  of  rotor  exit 
data  using  the  two  methods  are  compared.  In  order  to  verify 
the  experimental  procedures  at  low  cost^  a  bench-top  periodic 
flow  generate  :  was  used  first  to  provide  signals  similar  in 
type  and  frequency  to  those  from  the  Transonic  Compressor  in 
which  the  rotor  exit  flow  was  measured.  The  data  measured  in 
the  compressor  were  analyzed,  and  problems  encountered  in  the 
routine  application  of  the  reduction  methods  were  identified. 

The  possible  application  of  similar  techniques  in  the  Ro¬ 
tating  Annular  Cascade  (RGG)  at  the  DFVLR  Gottingen  is  as¬ 
sessed  and  recommendations  are  made  for  the  coordination  of 


future  experimental  programs . 


Accession  For 

WXS  GRAM 
DTIC  TAB 
Unannounced 
Justification- 


V 

V  V  V 


By - - - 

Distribution/ _ 

Availability  Codes 
Avail  and/or 
Dist  Special 


TABLE  OF  CONTENTS 


I.  INTRODUCTION .  10 

II.  REVIEW  OF  PUBLICATIONS  DEALING  WITH  UNSTEADY  FLOW 

MEASUREMENTS  IN  TURBOMACHINES  .  14 

III.  REVIEW  OF  THE  EXPERIMENTAL  EQUIPMENT .  16 

A.  TRANSONIC  COMPRESSOR  (TX)  .  16 

B.  INSTRUMENTATION .  17 

1.  Probes  for  Radial  Surveys .  17 

2.  Wall  Pressure  Measurement .  20 

3.  The  Timing  Disk .  21 

C.  DEVICES  FOR  DATA  ACQUISITION .  22 

1.  Hewlett  Packard  HP-21MX  Computer  .  22 

2.  Pacer .  22 

3.  Hewlett  Packard  HP-5610A  A/D  Converter  .  .  25 

IV.  DUAL-PROBE  DIGITAL  SAMPLING  (DPDS)  TECHNIQUE  ...  27 

A.  METHOD  I— ANALYTICAL  DATA  REDUCTION .  29 

B.  METHOD  II— NUMERICAL  DATA  REDUCTION .  32 

C.  ON-LINE  CALIBRATION  (OLC)  .  34 

V.  SOME  RESULTS  FROM  HIGH  RESPONSE  PRESSURE  PROBES  .  38 

A.  PRELIMINARY  INVESTIGATIONS  .  38 

1.  Results  from  Periodic  Flow  Generator  (PFG)  38 

2.  Calibration  of  New  Type  A  and  Type  B 

Probes .  40 

a.  Results  for  the  Type  A  Probe  .....  42 

b.  Results  for  the  Type  B  Probe .  45 


4 


3.  Application  of  Reduction  Programs  to 

Calibration  Data  . 

a.  Results  for  Analytical  Data  Reduction 

b.  Results  for  Numerical  Data  Reduction 

B.  RESULTS  FROM  THE  TX  . 

1.  Results  for  A  and  B  Probes  . 

2.  Results  for  Wall  Static  Pressure 

Transducers  . 

VI.  PROPOSAL  FOR  A  FIRST  APPLICATION  IN  THE  RGG  .  . 

VII.  CONCLUSION  . 

LIST  OF  REFERENCES  . 

FIGURES  . 

DISTRIBUTION  LIST  . 


' 

S 

i 


$ 


i 


♦t 

0> 


•I 

n 

i 


LIST  OF  FIGURES 


1.  Transonic  Compressor  (TX)  .  92 

2.  Geometry  of  the  Combination  Probe .  93 

3.  Arrangement  of  Probes  in  the  TX .  95 

4.  Probe  Geometries  .  96 

5.  Configuration  of  Probes  and  Locations  of  Case-wall 

Transducers .  97 

6.  DPDS- technique .  98 

7.  Data  Acquisition  System .  99 

8.  Schematic  of  "Pacer" . 100 

9.  Synchronized  Sampling  Signals  .  101 

10.  Results  from  the  PFG;  pn-i  0  _  *  151  in  H20; 

Pref  *  80  in  H20;  f  *  4‘.  5  kHz;  rep  »  10 . 102 

11.  Results  from  the  PFG;  Ppienum  a  146  in  H20; 

Pre£  *  80  in  HjO;  f  *  9.8  kHz;  rep  *10 . 103 

12.  Frequency,  generated  with  wave-generator;  f  *  10 

f  ■  10  kHz . 104 

13.  Calibration  facility  geometry  (not  to  scale)  .  .  .  105 

14.  pA  *  f  (a)  106 

15.  pA  »  f(a) . 107 

16.  pA  »  f  (a)  108 

17.  cPa/gP0A  versus  a;  M  =  .6 . 109 

18  *  CPA/CP0A  *  f(a/ag) . 110 

19.  CpQA  *  f($);  experimental  results . Ill 

20.  cPqa  *  after  curve-fit . 112 

21.  a  »  f(1>);  experimental  results . 113 

s 


*s 

1 

i 


6 


22.  a  *  f(<j>);  after  curve-fit . 114 

23.  pB  ■  f(o)  . . 

24.  cpB  *  f  (a)  ;  M  =  .6 . 116 

25‘  cP0B  *  117 

26.  Results  of  analytical  data  reduction  applied  to 

calibration  data  (old  A-  and  B-probes) . 118 

27.  Results  of  analytical  data  reduction  applied  to 

calibration  data  (old  A-  and  B-probes) . 119 

28.  Results  of  analytical  data  reduction  applied  to 

calibration  data  (new  A-  and  B-probes) . 120 

29.  Results  of  numerical  data  reduction  applied  to 

calibration  data . .  .  121 

30.  Results  of  numerical  data  reduction  applied  to 

calibration  data . 122 

31.  Results  for  the  A-  and  B-probe  with  different 

repetitions . 123 

32.  Probe  signals  taken  from  oscilloscope  .  124 

33.  Results  for  B-probe;  Rep  =  1 . 125 

34.  Results  for  B-probe;  Rep  *  10 . 1 26 

35.  Qualitative  explanation  of  wave  behavior  .  127 

36.  A  and  B  probe  signals;  Rep  =*10 . 128 

37.  A  and  B  probe  signals;  Rep  =  10 . 129 

38.  A  and  B  probe  signals;  Rep  =  40 . 130 

39.  A  and  B  probe  signals;  Rep  =40 . 131 

40.  pA  =  f(a)  for  A  probe  (step  1  to  128);  Rep  =  10  .  132 


41.  pA  =  f(a)  for  A  probe  (step  129  to  256);  Rep  =  10  133 

42.  pA  =  f(a)  for  A  probe  (step  35  to  64);  Rep  =  10  .  134 

43.  pB  »  f(a)  for  B  probe  (step  1  to  128);  Rep  =  10  .  135 

44.  pB  =  f(a)  for  B  probe  (step  129  to  256);  Rep  =  10  136 


7 


PB  »  fCcU  for  B  probe  (step  11  to  41);  Rep  =  10 

pA  »  f(a)  for  A  probe  (step  1  to  128);  Rep  =  40 

Results  from  A  probe;  a  =  52°  . 

Results  from  B  probe;  a  =  103°  . 

PA  »  f(oO  for  A  probe  (count  31) . 

M- number  and  flow  angle  distribution  versus 
peripheral  position  . 

Correction  for  a.  values  (position:  outside  wake 
flow) . . 

Correction  for  a,  values  (position:  within  wake 
flow)  . A . 

Flow  quantities  for  corrected  values  .  . 

Case-wall  results;  Rep  a  20  . 

Signals  from  K9  for  different  blade-pairs  .  .  . 

Curve-fit  for  signals  from  different  numbers  of 
repetitions  . 

Test  section  RGG  (a  =  45°)  . 

Proposal  of  trigger  system  for  RGG  . 

Sketch  of  total  pressure  probe  for  RGG  (not  to 
scale)  . 

Sketch  of  static  pressure  probe  for  RGG  (not  to 
scale)  ...  . 


ACKNOWLEDGMENT 


The  work  described  in  this  report  was  carried  out  while 
the  author  was  supported  for  twelve  months  by  the  DFVLR  and 
by  the  NPS.  The  work  was  completed  while  the  author  was  Naval 
Air  Systems  Command  Visiting  Research  Professor  in  Aeronautics 
from  October  1981  to  March  1982.  The  author  acknowledges  t 
financial  support  of  those  different  institutions. 

Furthermore,  the  author  would  like  to  express  his  many 
thanks  to  the  director  of  the  TPL,  Dr.  Raymond  P.  Shreeve, 
for  his  help  and  all  the  time  he  invested  in  the  technical 
discussions.  In  addition,  the  support  from  all  members  of 
the  TPL  is  also  appreciated  very  much.  Personal  friendship, 
human  integrity  and  technical  knowledge  of  the  people  involved 
in  the  project  made  the  l*a  year  period  an  unforgettable  time 
for  the  author. 


I .  INTRODUCTION 


It  is  one  important  goal  in  today's  axial  turbine  de¬ 
sign — for  stationary  steam  and  gas  turbines  in  electricity 
generating  plants  as  well  as  for  jet  engines— to  develop  ma¬ 
chines  which  are  as  compact  as  possible.  Therefore  larger 
enthalpy  differences  per  stage  are  used  in  today's  designs 
than  was  the  case  in  the  past.  This  results  in  supercritical 
pressure  ratios  for  the  conditions  occurring  in  some  parts  of 
the  blade  rows,  and  leads  to  flows  which  are  partly  subsonic 
and  partly  supersonic.  This  regime  is  usually  termed  tran¬ 
sonic,  and  is  well  known  to  be  difficult  to  investigate  both 
theoretically  and  experimentally.  This  holds  true  for  exter¬ 
nal  flows — flow  past  single  airfoils,  for  example,  but  it  is 
particularly  true  for  internal  flows — flow  within  turbomachines 
being  an  extreme  case. 

For  large  power-generating  steam  turbine  plants  the  shaft 
speed  is  fixed  by  the  RPM  of  the  generator.  Therefore  an  in¬ 
crease  in  work  output  is  only  possible  by  increasing  the  mass 
flow  rate.  As  a  result,  progressively  larger  diameters  of 
these  kinds  of  machines  are  being  used.  Near  the  tip  sections 
in  the  last  stages  of  low  pressure  steam  turbines  this  also 
results  in  transonic  relative  velocities. 

Presently,  it  is  not  possible  to  fully  describe  the  pas¬ 
sage  of  the  flow  through  the  different  stages  of  a  turbine 


from  either  experimental  or  theoretical  approaches.  In 
reality,  very  complicated  three  dimensional  unsteady 
flow  conditions  exist. 

However,  in  order  to  simplify  the  picture  for  design  pur¬ 
poses,  two  dimensional  flow  models  are  usually  used.  Experi¬ 
mentally,  two  dimensional  cascades  of  axial  blade  row  segments 
are  generated  either  in  an  annular  cascade  of  high  hub-to-tip 
ratio  or  by  developing  the  coaxial  stream  surface  of  interest 
into  a  plane  surface  as  a  rectilinear  cascade.  Because  of 
the  limited  number  of  blades  which  must  be  used,  the  simula¬ 
tion  of  the  required  periodic  flow  field  in  a  rectilinear 
cascade  is  usually  imperfect,  and  this  is  particularly  true 
for  transonic  flow  conditions. 

This  disadvantage  is  avoided  in  the  case  of  an  annular 
cascade,  and  for  some  years  such  a  facility  for  Rotating 
Annular  Cascades  (RGG)  has  been  in  operation  in  the  DFVLR- 
AVA  in  Gottingen,  West  Germany.  In  developing  this  facility, 
it  was  necessary  to  develop  a  new  method  of  evaluating  the 
measurements  such  that  the  results  were  representative  of  the 
flow  through  a  truly  two  dimensional  cascade.  The  method  was 
required  to  be  simple  but  reliable.  Such  a  method  was  de¬ 
veloped  and  was  verified  by  investigations  on  two  cascades  of 
very  different  types  [1,  2], 

The  fundamental  advantage  of  the  method  is  that  the  ex¬ 
perimental  determination  of  only  the  total  pressure  and 
total  temperature  in  the  homogenous  downstream  flow  field  is 


necessary.  This  requires  the  existence  in  the  downstream 
flow  field  of  a  region  where  nor-uniformities  from  wakes  have 
disappeared  and  where  the  boundary  layers  from  hub  and  tip  of 
the  annulus  have  not  penetrated.  Although  the  existence  of 
such  a  region  was  shown  it  was  considered  desirable  to  mea¬ 
sure  peripheral  and  time-varying  flow  properties  as  well. 

In  order  to  measure  the  unsteady  flow  field  behind  a 
rotor  one  must  either  use  high  response  probes  or  an  optical 
technique  of  some  kind.  While  laser  velocimetry  techniques 
have  many  advantages,  they  are  normally  able  to  measure  only 
two  components  of  the  velocity.  They  do  not  provide  measure¬ 
ments  of  the  pressure  field,  and  in  highly  unsteady  flows  the 
accuracy  must  be  questioned. 

At  the  Turbopropulsion  Laboratory  (TPL)  of  the  Naval 
Postgraduate  School  (NPS) ,  the  dual-probe  digital  sampling 
(DPDS )  technique  is  being  developed.  The  technique  is  based 
on  the  use  of  two  Kulite  semi-conductor  pressure  probes  of 
very  simple  design.  The  approach  has  the  potential  of  evalu¬ 
ating  all  three  components  of  the  periodic  velocity  vector  in 
a  turbomachine  or  rotating  cascade.  Following  early  work 
with  a  single  probe  in  a  compressor  flow  field  [21],  effort 
at  NPS  was  devoted  to  developing  and  verifying  a  method  of 
calibrating  and  applying  the  dual-probe  system  in  the  steady 
flow  of  a  free  jet  [3],  The  application  of  the  probe  system 
in  the  compressor  flow  field  was  delayed  until  recently  by 
drive  turbine  bearing  failure. 
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Now,  some  first  data  from  the  compressor  and  early  at¬ 
tempts  to  analyze  the  data  are  reported  here.  While  there  are 
some  open  questions  left  by  the  work  done  so  far,  it  is  recom¬ 
mended  nevertheless  that  the  method  be  applied  in  the  RGG  at 
the  DFVLR.  Because  of  the  early  stage  of  development  and 
complexity  involved  in  the  process  it  is  suggested  that  the 
efforts  at  the  DFVLR  and  NFS  be  conducted  in  close  association 
so  that  improvements  sure  directly  available  to  both  institutions. 

This  account  should  be  viewed  as  a  documentation  of  the 
learning  process  that  the  author  has  undergone.  It  reviews 
the  earlier  work  listed  as  references.  Key  passages  from 
those  references  have  been  included  without  change  on  occa¬ 
sion  when  any  contraction  of  the  wording  would  have  detracted 
from  the  descriptions  originally  given.  The  list  of  references 
is  organized  so  as  to  group  the  materials  logically,  rather 
than  in  the  order  that  they  are  mentioned  in  the  text.  The 
references  are  reviewed  in  the  following  section.  In  sections 
III  and  IV,  descriptions  of  both  the  experimental  and  computer 
hardware  are  given.  Section  V  describes  the  DPDS-method  and 
in  section  VI  some  results  are  shown.  Finally,  in  section  VII 
ideas  for  a  first  application  in  the  RGG  are  described. 


II.  REVIEW  OF  PUBLICATIONS  DEALING  WITH 


UNSTEADY  FLOW  MEASUREMENTS  IN  TURBOMACHINES 

In  [3]  to  [43]  publications  are  listed  which  are  a  useful 
source  of  information  in  the  field  of  unsteady  flow  measure¬ 
ments  in  turbomachines.  The  quoted  references  are  those  that 
the  author  has  used;  the  list  does  not  cover  the  available 
literature  completely.  It  should  be  mentioned  that  in  the 
following  text,  the  references  referred  to  most  often  are 
those  which  concern  the  research  activities  of  the  TPL. 

The  sequence  of  the  reference  list  is  such  that  Masters 
Theses  of  the  NPS  ([4]  to  [10])  are  given  first,  followed  by 
Technical  Reports  of  the  NPS  ([12]  to  [17])  and  finally, 
publications  in  the  open  literature  ([18]  to  [21]  and  [3]), 
which  generally  summarize  much  of  the  useful  material  from 
the  theses  and  reports.  In  [11]  a  U.S.  patent  is  referred 
to.  References  from  [22]  to  [43]  deal  mostly  with  different 
kinds  of  high  response  probes  used  in  the  flow  field  down¬ 
stream  of  rotors.  Some  results  from  laser  measurements  are 
reported  in  [37]  and  the  use  of  Schlieren  pictures  for  in¬ 
terpreting  the  flow  field  in  [40],  In  [31],  [39],  [41]  and 
[42]  the  use  of  high  response  pressure  transducers  to  measure 
the  dynamic  flow  behavior  on  rotor  blades  and  stator  blades 
are  reported.  Case  wall  measurements  based  on  the  same  kind 
of  instrumentation  are  reported  in  [36]  and  in  [32]  a  FORTRAN 
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program  is  given  for  reducing  case  wall  static  pressure  data 
to  contour  maps.  In  [24],  a  general  description  of  different 
problems  which  can  occur  in  investigating  flow  fields  in 
transonic  compressors,  is  given. 

The  ability  to  measure  the  true  time-averaged  static  and 
impact  pressures  is  of  considerable  importance  and  is  needed 
to  permit  a  full  interpretation  of  semi-conductor  probe 
measurements  in  non-steady  flow  fields.  This  is  discussed 
in  [26],  [27],  [33]  and  [34] .  For  application  in  transonic 
flow  fields  the  oil-filled  total  pressure  probe  described 
in  [33]  and  [34]  may  be  the  most  useful  tool  available. 


REVIEW  OF  THE  EXPERIMENTAL  EQUIPMENT 


A.  TRANSONIC  COMPRESSOR  (TX) 

The  transonic  compressor  is  located  in  the  second  test 
cell  of  the  High  Speed  Turbine  and  Compressor  Laboratory. 

This  laboratory  in  total  consists  of  three  test  cells,  a 
control  and  instrumentation  room  and  a  compressor  room  con¬ 
taining  a  single  12-stage  axial  air  compressor  that  provides 
the  power  for  all  experiments.  '  There  is  also  a  "spin  pit" 
for  the  determination  of  stresses  in  rotors.  Details  of  the 
laboratory  are  reported  in  [44]. 

Compressed  air  at  up  to  11  lbs.  per  sec.  at  three  atmos¬ 
pheres  pressure  can  be  valved  to  each  cell  from  an  Allis- 
Chalmers  multi-stage  axial  compressor.  The  1250  HP  compressor 
has  a  variable  speed  drive  with  an  automatic  control  and  a 
by-pass  surge  suppression  device.  The  mass  flow  rate  from 
the  compressor  is  metered  by  a  sharp  edged  orifice  meter  lo¬ 
cated  downstream  of  an  aftercooler  in  the  outlet  pipe. 

A  line  drawing  of  the  TX  is  shown  in  Fig.  1.  The  compres¬ 
sor  is  a  single  stage  axial  fan  designed  to  operate  at  30,460 
RPM  with  a  relative  tip  Mach  number  of  1.5.  The  rotor  diameter 
is  27.94  cms  (11  inches),  and  the  hub-to-tip  ratio  at  the 
rotor  face  is  0.5.  At  the  design  RPM  and  referred  flow  rate 
of  8.62  kg/sec  (19  lb/sec),  the  relative  flow  angle  at  the 
rotor  tip  is  65°  and  the  design  pressure  ratio  is  1.6.  The 
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compressor  rig  is  open  loop  with  an  air  filter,  throttle  and 
flow  nozzle  in  the  induction  pipe  and  with  the  compressor 
exhaust  going  to  atmosphere  via  a  radial  diffuser.  The  flow 
rate  is  controlled  by  an  electro-hydraulic  rotating  throttle 
plate  located  in  the  inlet  duct.  The  inlet  of  the  transonic 
compressor  is  designed  to  contract  and  to  rapidly  accelerate 
the  flow  in  order  to  generate  near  uniform  conditions,  with 
thin  boundary  layers,  at  the  rotor  face. 

The  compressor  test  rig  has  the  advantage  of  providing 
redundancy  in  measurements.  The  operating  power  can  be  de¬ 
termined  in  three  ways:  either  from  the  mass  flow  rate  and 
temperature  drop  through  the  turbine,  or  from  the  mass  flow 
rate  and  temperature  increase  through  the  compressor,  or  from 
the  speed  and  the  torque  measured  by  strain  gauges  on  flex¬ 
ures  between  the  compressor  stator  (mounted  on  bearings)  and 
the  machine  casing. 

B.  INSTRUMENTATION 

1.  Probes  for  Radial  Surveys 

Probe  survey  ports  are  provided  at  four  axial  planes 
and  fixed  Kiel  and  temperature  probes  are  installed  for  over¬ 
all  stage  performance  measurements  as  shown  in  Fig.  1.  The 
case  of  the  compressor  can  be  rotated  manually  about  the  axis 
while  the  machine  is  running  in  order  to  position  probes  and 
sensors  mounted  in  the  case  wall  so  as  to  obtain  surveys  in 
the  peripheral  direction  (with  respect  to  the  fixed  stator 
blading) . 


The  combination  temperature-pneumatic  probe  which  is 


used  is  shown  in  Fig.  2.  One  of  those  probes  is  used  in  front 
of  the  rotor  and  the  other  one  is  located  between  the  rotor 
and  the  stator,  in  use,  the  probes  must  be  rotated  about 
their  axis  to  balance  the  pressures  at  the  side  sensors,  p2 
and  p3,  and  the  yaw  angle  can  be  read  from  a  vernier  scale. 

The  center  sensor,  p^,  measures  an  indicated  impact  pressure, 
the  pressures  p2  and  p3  approach  static  pressure  and  the 
fourth  sensor,  p4,  provides  a  measurement  of  the  pitch  angle, 
when  the  probe  has  been  suitably  calibrated. 

The  probe  was  designed  specifically  to  be  used  down¬ 
stream  of  the  transonic  rotor  [4].  The  use  of  separate  sen¬ 
sors  (the  "claw"  arrangement  rather  than  pressure  taps  on  a 
shaped  surface)  was  to  minimize  the  "immersion"  effect.  It 
is  also  likely  that  the  separate  sensors  would  be  less  in¬ 
fluenced  by  unsteady  effects  from  the  rotor.  The  separation 
of  the  sensors  peripherally  is  of  no  consequence  because  sen¬ 
sors  at  the  same  radius  would  automatically  average  in  the 
same  passing  flow  field.  The  pitch  angles  in  the  compressor 
are  expected  to  be  small  and  severe  velocity  gradients  in  the 
radial  direction  are  not  expected  away  from  the  hub  and  tip 
regions.  The  necessary  radial  separation  of  the  fourth  sen¬ 
sor  is  therefore  not  considered  to  be  critical.  Four  sensors 
are  used  rather  than  five,  to  allow  the  simple  arrangement  of 
tubes  shown  in  Fig.  2.  This  arrangement  minimized  blockage 
in  a  strong,  simple,  leak-free  construction  and  gave  suitably 
short  response  times. 
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The  pneumatic  pressures  are  sensed  by  a  conventional 
differential  strain  gauge  transducer  and  Scanivalve  referenced 
to  atmosphere. 

The  separation  of  the  thermocouple  sensor  from  the  im¬ 
pact  pressure  sensor  in  the  radial  direction,  Fig.  2,  is  also 
not  important  since  it  could  be  properly  accounted  for  in  the 
data  reduction  from  radial  surveys. 

Probe  calibration  procedure  and  results,  data  from  ra¬ 
dial  surveys  and  stage  performance  measurements  are  reported 
in  [17]. 

In  [21]  a  method  was  described  for  using  two  semi¬ 
conductor  pressure  probes  together  with  the  technique  of  syn¬ 
chronized  sampling,  to  obtain  the  distribution  of  the  velocity 
vector  downstream  of  a  rotor.  The  geometries  of  the  two  probes, 
designated  Type  A  and  Type  B,  and  their  installation  in  the 
compressor  annulus  are  shown  in  Fig.  3.  The  two  probes  are 
positioned  as  shown  at  the  same  radial  displacement,  but  sepa¬ 
rated  peripherally  around  the  machine  annulus.  The  probes  are 
mounted  in  such  a  way  that  they  can  be  rotated  about  an  axis 
through  the  probe  tip.  The  details  of  the  geometry  of  the  two 
probes  are  shown  in  Fig.  4. 

The  particular  shape  of  the  B  probe  is  chosen  to  enable 
measurements  to  be  made  close  to  the  case  wall,  and  therefore 
limitations  occur  in  the  measurements  which  can  be  made  close 
to  the  hub. 
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It  was  argued  that,  in  principle,  by  rotating  the 
probes  in  yaw  about  their  tips  and  controlling  the  sampling 
of  the  data  from  each  probe  to  be  at  the  same  position  in  the 
rotor  frame,  the  system  of  two  separate  probes  could  be  used 
to  acquire  data  at  a  single  point  in  the  periodic  flow  from 
the  rotor,  thus  obtaining  data  corresponding  to  those  obtained 
from  the  sensors  of  four-  or  five-hole  pneumatic  probes  when 
measuring  velocity  in  steady  flows.  The  technique  involved 
the  use  of  probes  having  the  simplest  geometry  and  thus  avoided 
the  large  size,  expense  and  unreliability  of  multiple  sensor 
probes  which  incorporate  multiple  semiconductor  transducers. 

It  would  be  far  cheaper  to  replace  a  single  probe  of  simple 
geometry  in  the  event  of  a  failure  compared  to  repairing  or 
replacing  a  multiple  sensor  probe.  Because  of  the  simple  sen¬ 
sor  tip  geometry  (that  of  a  cylinder  at  incidence  to  the  flow) 
the  unsteady  response  was  likely  to  be  as  good  as  could  be  ex¬ 
pected  of  any  single  physical  sensor. 

All  probes  are  yawed  and  moved  to  different  radial 
positions,  manually.  Both  yaw  angle  and  radial  position  are 
recorded  automatically  using  position  potentiometer  circuits. 

2.  Wall  Pressure  Measurement 

Twelve  Kulite  CQL-080-25  pressure  transducers  with 
natural  frequency  about  125  kHz  are  mounted  with  their  dia¬ 
phragms  flush  with  the  inner  case  wall  of  the  compressor. 

The  transducers  are  located  along  four  axial  lines  separated 
peripherally  by  10  and  90  degree  arcs  as  shown  in  Fig.  5. 
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The  transducers  are  used  in  conjunction  with  Datel  Model  201C 
instrumentation  amplifiers  which  have  a  flat  frequency  response 
to  100  kHz.  Each  Kulite  pressure  transducer  is  matched  by  a 
pneumatic  static  pressure  tap  at  the  same  axial  location  in 
the  case  wall,  but  displaced  circumferentially.  Other  pneu¬ 
matic  static  pressure  taps  are  provided  upstream  and  downstream. 
The  reference  tubes  of  the  transducers  are  connected  by  separate 
plastic  tubes  to  a  manifold  fed  from  a  controlled  pressure 
source  or  to  the  pneumatic  line  from  the  pressure  tap  at  the 
same  axial  location.  Clamps  were  used  to  valve  the  reference 
side  of  the  transducer  either  to  the  regulated  pressure  or  to 
the  corresponding  pneumatic  pressure  tap. 

3.  The  Timing  Disk 

To  enable  synchronization  of  the  sampling  of  the  semi¬ 
conductor  pressure  transducer  outputs  with  the  rotation  of  the 
rotor,  a  timing  disk  is  fitted  to  the  rotor  shaft.  Fig.  6.  The 
disk  has  a  circle  of  eighteen  equally  spaced  holes  (equal  to 
the  number  of  rotor  blades)  at  one  radius  and  a  single  hole  at 
a  smaller  radius.  Light  emitting  diodes  and  photo  detectors 
aligned  with  the  holes  generate  bell-shaped  electrical  pulses 
at  one-per-revolution  and  one-per-blade  frequencies.  The  sys¬ 
tem  is  described  in  detail  in  [6]. 

The  blade  passing  frequency,  f,  can  be  computed  by 
multiplying  the  wheel  speed  (revolutions  per  second)  with  the 
number  of  blades.  The  maximum  rpm  of  the  TX  is  about  500/sec 
and  the  rotor  constsis  of  18  blades.  Therefore,  the  expected 
frequency  range  is  up  to  10  kHz. 
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C.  DEVICES  FOR  DATA  ACQUISITION 

The  primary  hardware  devices  used  for  data  acquisition  are 
the  HP-21MX  computer  with  its  20  Megabyte  magnetic  disc,  print 
er,  plotter,  and  two  terminals,  the  "Pacer"  and  the  HP-5610A 
16-channel  A/D  converter.  The  complete  data  acquisition  sys¬ 
tem  is  shown  in  Fig.  7. 

1.  Hewlett  Packard  HP-21MX  Computer 

The  HP-21MX  is  a  (Micro-programmable)  mini-computer 
having  128  machine  instructions  and  32K  of  logical  in  24QK  of 
main  frame  memory.  A  20  megabyte  capacity  disc  and  disc  drive 
surmount  the  core  memory  limitation  and  provide  multiple  use 
under  time  sharing. 

An  important  feature  which  is  typical  of  computers  of 
this  size  is  the  input-output  structure.  With  a  limited  num¬ 
ber  of  relatively  slow  I/O  devices  to  be  serviced,  the  com¬ 
puter  can  communicate  with  all  devices  through  a  single  port 
known  as  the  I/O  bus.  Each  device  requires  its  own  I/O  inter¬ 
face  on  the  bus.  The  interface  acts  as  a  filter  and  ensures 
that  output  information  is  received  only  by  the  device  desig¬ 
nated  to  receive  it  and  that  input  information  is  put  on  the 
bus  from  only  one  device  at  a  time.  The  I/O  software  drivers 
control  the  I/O  hardware  interfaces  by  commands  to  either 
"turn-on"  or  "turn-off". 

2.  Pacer 

This  peripheral  device  which  was  designed  and  con¬ 
structed  to  trigger  data  acquisition  from  a  stationary  sensor 
at  any  fixed  point  in  the  rotating  rotor  frame,  independent 
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of  the  rotor  speed,  is  shown  schematically  in  Fig.  8.  In  its 
original  form,  a  detailed  description  of  the  internal  opera¬ 
tion  is  given  in  [ 6 ] . 

The  bell-shaped  optical  signals  from  the  timing  disk 
are  first  converted  in  the  "wave  shaper"  section  into  square 
pulses  suitable  for  TTL  (Transisitor-Transistor-Logic)  circuits 

The  "Pacer"  consists  of  two  major  sections,  a  "synchro¬ 
nized  command  pulse  section”  and  a  "RPM  counting  section". 

The  "synchronized  command  pulse"  section  is  the  heart 
of  the  "Pacer".  The  one-per-rev.  signal  is  used  as  a  master 
pulse  in  the  timing  process.  The  one-per-blade  signal  is 
halved  in  frequency  through  a  counter  and  input  into  a  PLL 
(Phase-Locked  Loop)  circuit.  The  PLL  "locks-on"  and  outputs 
a  s imi lar  wave  form  at  256  times  the  input  frequency.  A  com¬ 
parator  compares  the  number  of  counts  following  the  master  re¬ 
set  (taken  from  the  one-per-rev.  signal)  with  the  programmed 
number  received  from  the  computer.  When  they  agree,  the  com¬ 
mand  signal  from  the  computer  to  the  A/D  converter  is  allowed 
to  pass  and  conversion  occurs. 

The  device  can  be  programmed  to  time  the  data  conver¬ 
sion  to  occur  at  any  point  in  the  rotor  frame,  since  the  "lock- 
on"  makes  it  independent  of  the  rotor  speed. 

The  pacer  can  also  be  programmed  to  allow  conversion 
either  at  a  single  location  in  a  particular  blade  passage,  or 
at  the  corresponding  positions  in  each  blade  passage. 

The  "RPM  counting  section"  continuously  counts  the  num¬ 
ber  of  250  kHz  time  base  pulses  (obtained  from  the  computer's 
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1  MHz  clock)  that  occur  between  the  one-per -revolution  pulses 
received  from  the  test  rig.  This  number  of  counts  is  available 
as  an  output  from  the  '’Pacer"  on  every  revolution  of  the  rotor. 

The  original  "Pacer"  was  recently  modified  to  overcome 
the  following  three  limitations  which  were  realized  when  using 
the  original  version: 

(1)  The  timing  pulses  generated  with  the  "Pacer"  were 
not  always  spaced  linearly  in  time  between  blade  pair  synchro¬ 
nizing  pulses.  This  resulted  in  data  which  in  some  cases  was 
subtly  distorted,  and  in  other  cases  appeared  to  have  noise 
riding  on  it. 

(2)  The  range  over  which  the  "Pacer"  could  follow 
rotor  RPM  changes  and  remain  synchronized  was  limited  to  ap¬ 
proximately  ±15%  of  the  initial  RPM  at  which  the  "Pacer"  was 
set  to  take  the  data.  This  required  a  careful  manual  adjust¬ 
ment  of  trimming  capacitors  and  resistors  to  obtain  "lock-on" 
in  each  range  of  RPM  in  order  to  acquire  data  at  different 
speeds. 

(3)  The  rate  at  which  data  could  be  taken  was  limited. 
Rather  than  being  able  to  sample  data  on  every  revolution  of 
the  rotor,  the  system  was  only  capable  of  taking  data  once 
every  8  to  10  revolutions,  depending  on  RPM. 

Changes  to  the  "Pacer"  are  described  in  detail  in  [10] 
which  fully  overcame  the  above  three  limitations.  They  in¬ 
volved  changes  in  both  hardware  and  software  over  those  re¬ 
ported  in  [6].  The  present  system  is  extremely  easy  to  use 
and  presents  no  obvious  limitations. 


.  Hewlett  Packard  HP-5610A  A/D  Converter 


Up  to  sixteen  channels  of  analog  input  voltages  in  the 
range  of  ±1.0  volt  can  be  sampled  and  recorded  by  the  A/D  con¬ 
verter  under  the  control  of  the  21MX  computer.  The  acquisition 
process  is  governed  by  the  built-in  features  of  the  A/D  con¬ 
verter  and  by  the  programs  which  can  be  written  for  the  21MX. 

The  HP-5610A  Analog-to-Digital  Converter  is  designed 
to  convert  analog  voltages  into  ten  data  BIT'S  at  the  maximum 
throughput  rate  of  100  kHz.  The  HP-5610A  can  operate  in  one 
of  six  modes.  The  two  fundamental  modes  of  operation  are  "se¬ 
quential" — in  which  all  16  channels  are  scanned  in  numerical 
order,  and  "random  access" ,  in  which  a  specific  channel  is 
sampled  on  receipt  of  a  command  word  and  an  encode  command 
pulse  from  the  21MX  computer. 

Variations  within  these  two  basic  modes  are  also  pos¬ 
sible  which  allow  for  externally  supplied  sequencing  pulses 
and  for  externally  controlled  timing  of  samples. 

Another  mode  which  is  used  only  for  non-paced  data  is 
the  "Free  Run,  Random  Access"  mode.  In  this  mode  the  command 
word  is  required  as  before,  but  no  encode  command  is  needed 
from  the  computer.  The  A/D  converter  simply  converts  data  as 
fast  as  it  can  (100,000  samples  per  second)  on  the  selected 
channel. 

The  command  word  tells  the  A/D  converter  which  mode  of 
operation  to  use  and  which  channel  number  to  sample.  The  en¬ 
code  command  pulse  triggers  the  data  conversion  to  start  2  usee 
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later.  The  data  conversion  itself  is  finished  in  a  total  time 
of  10  usee.  Using  computer-issued  encodes ,  which  is  the  mode 
required  for  paced  data,  the  sample  cycle  time  is  2Q  usee. 
Hence  data  can  be  converted  at  rates  of  up  to  50,000  samples 
per  second,  depending  on  how  rapidly  each  successive  command 
word  is  received.  With  an  input  conversion  aperture  of  50 
nanoseconds,  rapidly  changing  signals  (100  kHzl  can  be  con¬ 
verted  accurately. 

The  A/D  digitizes  data  using  a  method  of  successive 
approximations.  That  is,  a  known  test  voltage  is  compared 
with  the  sample  voltage  beginning  with  1/2  full  scale  (0.5v 
in  the  present  system)  and  a  decision  is  made  whether  the  test 
voltage  is  high  or  low.  The  comparisons  continue,  halving  the 
range  each  time  until  the  machine  accuracy  is  reached.  For  a 
resolution  of  ±1.0  volt  to  10  BIT'S  the  accuracy  of  the  A/D 
converter  is  1/512  =  1.95  mV. 


IV.  DUAL-PROBE  DIGITAL  SAMPLING  (DPPS)  TECHNIQUE 


Until  the  development  of  miniature  semiconductor  trans¬ 
ducers  having  large  bandwidth  frequency  response#  real-time 
pressure  measurements  in  high  speed  compressors  were  not 
practicable.  Currently,  using  these  pressure  transducers, 
coupled  with  fast  analog-to-digital  conversion  and  computer 
controlled  time  delay  circuits,  a  variety  of  new  real-time 
measurement  techniques  are  possible.  The  well-known  tempera¬ 
ture  dependence  of  the  transducers  (see  for  example  [39]), 
can  be  overcome  by  an  on-line  calibration  (OLC)  such  as  is 
explained  in  Section  IV  C. 

The  difficulties  to  be  overcome  in  the  design  of  a  sam¬ 
pling  method  include  achieving  timing  accuracy  at  high  blade 
passing  frequencies  and  accounting  for  the  small  variations 
in  rotor  speed  that  are  always  present. 

The  technique  of  DPDS  applied  in  the  transonic  compressor 
is  illustrated  in  Fig.  6  and  Fig.  9.  Referring  to  these 
figures,  a  transducer  which  is  positioned  downstream  of  the 
rotor  and  fixed  to  the  machine  frame  outputs  a  time  varying 
signal  which  is  basically  periodic  at  the  blade-passing  fre¬ 
quency.  By  successively  triggering  the  sampling  of  the  trans¬ 
ducer  with  a  fixed  delay  following  a  one-per-rev.  trigger 
pulse,  the  condition  at  a  fixed  point  in  the  rotor  frame  can 
be  determined.  By  programming  the  delay  time,  in  principle 


any  chosen  location  in  the  rotor  frame  can  be  examined.  How¬ 
ever,  at  design  speed,  the  blade  passing  frequency  is  near  10 
kHz,  and  conventional  time  delay  techniques,  in  the  presence 
of  speed  fluctuations,  can  not  provide  the  required  resolution. 
A  computer  peripheral  device  called  "Pacer"  was  therefore  de¬ 
signed  and  developed  to  provide  the  programmable  delay  ([6] 
and  [10])  . 

If  the  stagnation  temperature  is  known,  four  quantities 
are  sufficient  to  determine  the  velocity  vector:  for  example, 
static  pressure,  pw,  total  pressure,  pfc,  yaw  angle,  a,  and 
pitch  angle,  $.  Following  calibration  however,  any  four  inde¬ 
pendent  measurements  of  pressures  can  be  used  to  evaluate  these 
same  unknowns  (as  is  done  typically  with  a  multiple  sensor 
probe  in  steady  flow  [17]).  In  principle,  the  required  mea¬ 
surements  can  be  obtained  by  rotating  the  probes  about  their 
tips  to  four  different  directions  where  pressure  readings  can 
be  taken.  In  reality,  one  must  pay  attention  to  the  character¬ 
istic  behavior  of  the  probes,  i.e.  the  output  of  the  probes  as 
a  function  of  yaw  angle,  pitch  angle  and  Mach  number  (M)  in 
the  ranges  of  these  parameters  which  are  to  be  measured  in  the 
compressor  flow  field.  Experience  shows  in  fact  that  it  is 
preferable  to  evaluate  curves  of  transducer  pressure  (p)  as  a 
function  of  probe  yaw  angle  setting  (01^)  and  to  refer  these  to 
calibration  curves.  It  is  obvious  that  probe  readings  at  only 
four  different  yaw  positions  can  not  provide  a  precisely  de¬ 
fined  curve  and  therefore  a  comparison  which  is  sufficiently 
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accurate.  Certainly ,  more  than  four  different  probe  settings 
sure  necessary.  Seven  positions  were  suggested  in  [3];  nine 
and  eleven,  respectively,  are  used  for  data  reported  here. 

In  principle,  there  are  many  possible  methods  for  reducing 
data  acquired  in  the  DPDS  technique.  To  date  three  entirely 
different  approaches  have  been  attempted  at  the  TPL.  The 
first,  Method  I,  uses  a  particular  set  of  polynomial  approxi¬ 
mations  for  the  characteristics  of  the  two  probes  obtained  by 
calibration  [3].  The  second,  Method  II,  is  described  in  [14], 
and  requires  the  representation  of  the  calibration  as  an  array 
of  pressure  coefficients.  The  method  is  a  further  development 
of  the  work  reported  in  [13]  and  overcomes  some  of  the  earlier 
limitations.  The  description  and  results  of  the  third  method. 
Method  III,  which  makes  use  of  work  reported  in  [15]  in  an  ap¬ 
proach  which  parallels  that  used  for  steady-state  velocity 
probes  [17],  is  given  in  [45], 

In  this  report,  TX  data  are  analyzed  using  Method  I  and 
Method  II  only.  Some  details  of  the  procedures  will  be  given. 

A.  METHOD  I— ANALYTICAL  DATA  REDUCTION  [3] 

Using  the  probe  pressure  coefficient,  (c  ) ,  defined  as 


where  p  *  pressure  at  the  probe  sensor,  p  =  static  pressure 
and  pfc  *  total  pressure,  the  characteristic  behavior  of  the 
A-probe  as  a  result  of  observation  (first)  and  calibration 


(second)  can  be  analytically  expressed  by  the  following  three 
fourth  order  polynomial  equations: 
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where  aA  is  the  flow  yaw  angle  with  respect  to  the  axis  of 
the  probe#  agA  indicates  the  value  of  aA  where  static  pres¬ 
sure  is  measured  and  $  is  the  pitch  angle.  cp0A  defines  the 

maximum  c  value . 

P 

For  the  B-probe,  only  the  equivalent  to  equation  (3)  is 
useful  because  its  inclination  (see  Fig.  4)  is  such  that  the 
angle  at  which  static  pressure  is  registered,  asB»  does  not 
exist  at  most  positive  angles  of  pitch.  The  fourth  expression 
is  therefore 
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A.^,  C^,  and  (i  =  0  to  4)  are  polynomial  coefficients 
which  are  established  by  calibration.  In  Section  V  some  more 
details  are  given  when  the  evaluation  of  these  curves  for  two 
new  A-  and  B-probes  is  described. 

In  the  above  equations  for  every  probe  yaw  angle  position 


which  is  set  in  the  data  acquisition  procedure  (up  to  11 


positions  have  been  used)  a  single  value  for  each  of  a,  <J>,  pg 
and  pfc  can  be  derived. 

The  flow  yaw  angle,  a,  is  taken  simply  as  the  position  of 
the  maximum  of  the  probe  output  when  the  probe  pressure  is 
plotted  versus  the  probe  yaw  angle  setting  in  the  machine, 
oipA(m) .  Of  course,  this  could  be  done  for  the  output  of  both 
probes  and,  in  the  ideal  case,  the  same  value  would  be  cal¬ 
culated.  To  apply  the  calibration  equations  the  yaw  angles 
must  first  be  calculated  relative  to  the  axis  of  the  probe; 
i.e.  aa(m)  *  oipA(m)  -  a. 

The  pitch  angle,  <t> ,  can  then  be  calculated  iteratively  for 
every  data  point,  m,  from  the  following  implicit  equation  for 
<J>  obtained  using  Eq.  (1)  and  Eq.  (2) : 
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In  Eq.  (6) 


(7) 


where  cpQA  and  cpo3  are  given  by  Eq.  (3)  and  Eq.  (5)  respec¬ 
tively,  and  agA  is  given  by  Eq.  (4).  pQA  and  pQB  are  the 

maximum  values  evaluated  by  fourth  order  curve  fits  through 
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the  m  data  points,  for  the  type  A  and  type  B  probes  respec¬ 
tively.  It  should  be  noted  that  only  the  maximum  value  of  the 
B-probe  is  used  in  this  procedure.  The  solution  of  Eq.  (6)  to 
obtain  can  be  carried  out  using  Newtonian  iteration. 

The  static  pressure  for  every  data  point,  p  (m) ,  can  be 
calculated  using 


p3(m) 


Pqa  -  Pqb'1|'u1 

1  -  tp(<P) 


(8) 


By  applying  Eq.  (1)  to  pQA  the  total  pressure  for  every  data 
point,  pt(m),  can  be  obtained  using 
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It  can  be  seen  from  Eq.  (6)  that  for  pA(m)  *  pQA  the  left 
side  of  the  equation  becomes  unity.  In  the  ideal  case  (per¬ 
fectly  symmetric  probe  tip  geometry)  ^(m)  should  be  equal  to 
0  when  pA(m)  »  pQA,  aq  *  and  no  solution  for  <Hm)  is  then 
possible. 


B.  METHOD  II — NUMERICAL  DATA  REDUCTION  [14] 

Method  II  uses  the  values  of  c  as  defined  by  Eq.  (1)  and 

p 

obtained  from  calibration  data  in  the  form  of  arrays  corre¬ 
sponding  to  a  particular  grid  of  yaw  and  pitch  angles.  Since 
stem  effects  and  manufacturing  imperfections  will  lead  to  de¬ 
partures  from  the  perfectly  symmetric  probe  behavior  in  yaw, 
a  method  which  does  not  require  an  analytical  expression  for 
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the  probe  characteristics  has  potential  advantages.  The 
authors  of  [14]  developed  a  numerical  solution  which  avoided 
any  requirement  for  curve  fitting  and,  in  principle,  allowed 
a  determination  of  ps,  pf, ,  a  and  <t>  with  only  four  measure¬ 
ments.  In  the  current  version,  the  pressure  coefficients 
defined  in  Eq.  (1)  are  assumed  to  be  totally  insensitive  to 
Mach  number  and  Reynolds  number.  Nevertheless,  if  any  sig¬ 
nificant  variations  did  occur  an  additional  iteration  scheme 
could  be  added  to  the  computer  solution. 

From  four  different  measurements— three  from  the  A-probe 
and  one  from  the  B-probe  are  suggested — four  c  values  can  be 

Jr 

calculated.  It  is  assumed  that  each  characteristic  has  only 
one  maximum  in  the  range  of  interest  and  therefore  one  com¬ 
bination  of  yaw  and  pitch  angle  (.a,  $)  exists  where  the  four 
Cp  values  fit  those  in  the  calibration  tables.  The  FORTRAN 
computer  program  given  in  [14]  provides  a  systematic  trial- 
and-error  search  process,  essentially  a  convergence  scheme 
on  the  two  variables  yaw  angle  and  pitch  angle.  A  linear, 
double- interpolation  scheme  is  available  as  well  as  a  second- 
order,  double-interpolation  scheme. 

The  computer  program  in  [14]  was  entered  into  the  HP-21MX 
computer  together  with  the  c  data  arrays,  also  given  in  [14], 

hr 

which  are  based  on  analytical  functions.  The  Cp  values  so 
derived  are  automatically  symmetric  for  positive  and  negative 
yaw  angles,  a. 


First,  the  program  was  run  for  the  examples  given  in  [14]. 
The  results  agreed  perfectly.  Second,  the  method  was  applied 
using  the  actual  (unsymmetric)  c^  calibration  data  arrays. 
Incomprehensible  and  obviously  incorrect  results  occurred. 
After  changing  row  and  column  indices  of  the  c  arrays  in  the 

r 

FORTRAN  program  (referring  to  [14],  statements  3010  to  3040, 
page  31  and  statements  4830  to  4900,  page  44;  in  addition, 
IFLAG  =■  0  is  omitted  after  statement  4990)  results  were 
evaluated,  which  were  in  the  expected  range.  On  rerunning 
the  corrected  program  using  the  example  data,  only  minor  de¬ 
viations  from  the  first  results  were  found.  Hence  the  cor¬ 
rected  program  is  thought  to  be  working  for  non-symmetric 
probe  characteristics.  Additional  checks  using  calibration 
data  from  two  new  probes  are  reported  in  Section  V. 

C.  ON-LINE  CALIBRATION  (OLC) 

The  reduction  procedures  described  in  the  previous  two 
sections  require  the  input  of  correct  values  of  the  probe 
pressures.  Because  of  the  temperature  sensitivity  of  the 
transducers,  which  affects  mainly  the  zero  balance  but  to  a 
lesser  degree  the  output  slope  of  the  strain  gauge,  correc¬ 
tions  have  to  be  made  to  the  transducer  voltages  to  obtain 
reliable  quantitative  pressure  measurements.  The  problem 
occurs  even  though  a  temperature  compensation  module  is  in¬ 
cluded  in  the  transducer  circuit.  The  problem  of  temperature 
dependence  is  widely  discussed  in  the  literature;  for  example, 
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[8],  [23],  [33],  [39].  Two  approaches  are  in  common  use; 
either  the  signals  are  used  as  a  measure  of  the  fluctuations 
only  and  the  time-averaged  pressure  level  is  measured  by  a 
second  device,  or  flow  temperature  is  measured  simultaneously 
and  corrections  are  computed  based  on  an  earlier  temperature 
calibration. 

Progressive  improvements  have  been  made  by  the  transducer 
manufacturers  to  reduce  temperature  sensitivity  during  the 
past  decade.  Also,  in  one  case,  automatic  electronic  com¬ 
pensation  similar  to  linearization  of  hot  wire  circuits,  is 
reported  [46]. 

In  the  TX  an  on-line  calibration  procedure  for  all  Kulite 
transducers  (A-  and  B-probe  and  wall  statics)  is  used  to  ac¬ 
count  for  temperature  dependent  effects.  For  the  A-probe  the 
method  is  described  in  detail  in  [8]  and,  therefore,  only  a 
brief  summary  is  given  here. 

To  evaluate  the  slope  (differential  pressure  vs.  voltage 
output) ,  x^,  of  the  transducers  in  situ,  while  ready  to  take 
data,  the  transducers  sure  referenced  to  different  pressures 
(6  at  maximum)  and  a  linear  least  squares  routine  is  used  to 
obtain  x^.  Experience  shows  that  almost  without  exception 
the  slope  is  linear  and  reproducible  to  within  the  accuracy 
of  measurement  for  the  modest  changes  in  temperature  involved 
in  the  experiments  (~10-20°C) . 

Additionally,  the  above  procedure  establishes  an  inter¬ 
cept  value,  Xg,  which  can  be  monitored  over  longer  times  to 
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examine  any  temperature  drift  which  may  occur  while  test  con¬ 
ditions  are  held  constant.  Experience  shows  that  normally 
deviations  equivalent  to  less  than  .5  in  H20  pressure  occur 
while  dynamic  pressures  are  never  less  than  20  in  H20.  Larger 
deviations  can  be  interpreted  as  poor  temperature  compensa¬ 
tion,  which  was  the  case,  on  one  occasion,  for  one  of  the  case 
wall  transducers  which  was  subsequently  replaced. 

In  order  to  interpret  the  intercept,  x^,  the  assumption 
is  made  that  a  pneumatic  probe  measures  the  correct  time 
averaged  pressure  when  positioned  at  the  same  radial  location 
and  set  to  the  same  yaw  angle  as  the  high  response  probe. 

This  is  to  be  regarded  as  a  delicate  problem  because  investi¬ 
gations  are  published,  for  example  [33],  [34],  which  clearly 
indicate  that  this  can  be  a  totally  wrong  assumption.  To 
date,  it  has  been  of  minor  importance  for  the  TX  measurements 
due  to  low  peak-to-peak  pressure  amplitudes  which  have  been 
measured  so  far.  An  example  is  given  in  [17]  and  there  it  is 
indicated  that  the  parameter  "ratio  of  wave  amplitude  to  mean 
value,"  on  which  the  error  depends,  is  10  times  less  than  the 
lowest  value  investigated  in  [34].  For  future  measurements 
at  higher  rotational  speeds,  care  will  be  taken  to  reexamine 
the  importance  of  averaging  errors  under  the  then  prevailing 
machine  conditions. 

When  the  transducers  are  referenced  to  the  time  averaged 
pneumatic  pressure  their  outputs  are  set  to  "0”.  In  applying 
the  calibration,  samples  of  Kulite  output  voltage  are  reduced 


to  values  of  pressure,  Pprobe#  using  the  calibration  coeff- 
cients  xQ  and  x,  in  the  equation 

Pprobe  ‘  x0  +  *1  '  e  +  Pre£  C10) 

where  xQ  is  expected  to  be  very  close  to  zero.  In  acquiring 
data,  the  transducer  reference  tube  is  usually  open  to  atmos¬ 
phere  and  the  reference  pressure,  pref,  therefore  equals  at¬ 
mospheric  pressure. 

The  total  pressure,  pfc,  read  by  the  combination  probe, 

Fig.  2,  is  used  for  the  A-probe  to  evaluate  the  intercept, 
xQ.  The  intercept  for  the  B-probe  can  not  be  found  in  this 
way  because  of  its  inclination.  Fig.  4.  Rather,  the  static 
pressure,  pg,  and  flow  pitch  angle,  <t>,  from  the  combination 
probe  are  used  in  Eq.  (5)  to  compute  the  maximum  output  of 
the  B-probe,  P0B.  This  procedure  is  necessary  to  adjust  the 
transducer  of  the  B-probe  to  "0"  on-line.  It  is  noted  that 
the  process  does  not  involve  equating  to  a  pressure  which  is 
measured  directly  in  the  TX.  In  principle  a  pneumatic  equi¬ 
valent  to  the  B-probe  geometry  could  be  used,  but  it  is  feared 
that  minor  differences  of  the  inclination  angle  or  tip  geome¬ 
try  could  lead  to  totally  unacceptable  errors.  The  present 
procedure  is  one  which  ensures  that  the  time*  ^raged  results 
of  the  Type  A  and  B  probe  system  are  the  same  as  those  given 
by  the  combination  probe.  Within  this  constraint,  the  Type  A 
and  B  system  can  resolve  the  time-dependent  behavior  as  a 
blade-to-blade  distribution  of  flow  properties. 
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V.  SOME  RESULTS  FROM  HIGH  RESPONSE  PRESSURE  PROBES 


A.  PRELIMINARY  INVESTIGATIONS 

1.  Results  from  Periodic  Flow  Generator  (PFG) 

A  small  periodic  flow  generator  (PFG)  was  used  by 
the  author  as  a  means  of  becoming  acquainted  with  the  high 
speed  data  acquisition  system  in  the  simplest  way.  The  ap¬ 
paratus  is  described  in  detail  in  [5]  and  [7]  and  only  the 
main  features  are  given  here. 

Regulated  air  is  delivered  through  a  settling  cham¬ 
ber  to  a  small  nozzle.  A  sprocket  (or  interrupter)  wheel 
driven  by  a  speed-regulated  electric  motor,  is  mounted  di¬ 
rectly  above  the  jet  such  that  the  teeth  of  the  sprocket  pass 
through  the  jet.  The  spaces  between  the  teeth  simulate  blade 
passages  and  "one  per  revolution"  and  "one  per  blade"  signals 
are  generated  from  a  similar  arrangement  of  holes,  light 
emitters  and  detectors  as  is  used  on  the  TX.  With  60  teeth 
on  the  wheel  and  10,000  rpm  from  the  motor,  blade  passing 
frequencies  up  to  10  kHz  are  obtained  to  simulate  the  condi¬ 
tions  expected  in  the  TX.  The  high  response  pressure  trans¬ 
ducer  is  mounted  on  the  centerline  of  the  jet  behind  the 
plane  of  the  sprocket  wheels,  so  that,  ignoring  transient 
flow  deflections  in  the  unsteady  process,  the  jet  impact 
pressure  is  measured.  In  the  potential  core  of  the  free 


jet,  the  transducer  would  be  expected  to  register  a  pressure 
equal  to  the  supply  pressure. 

In  addition  to  tests  conducted  solely  for  the  pur¬ 
pose  of  learning,  several  small  experiments  were  carried 
out  of  some  general  interest.  For  example,  results  were 
obtained  to  document  how  well  the  "Pacer"  works.  In  Figs . 

10  and  11  two  examples  are  shown  which  are  results  obtained 
to  verify  that  the  "Pacer”  was  performing  correctly  as  a 
trigger  device.  On  the  interrupter  wheel  one  "blade"  pair 
was  distinguished  by  fixing  a  wire  across  the  gap  between 
two  teeth.  By  using  the  "Pacer”  to  obtain  measurements  be¬ 
hind  the  different  "blade"  pairs,  the  blade  space  contain¬ 
ing  the  wire  would  be  detected  by  the  appearance  of  a 
"wake".  From  the  two  figures,  which  show  the  results  for 
two  different  frequencies,  it  can  be  seen  that  the  wake  was 
present  only  for  pair  2.  All  the  other  signals  were  quite 
similar,  as  also  were  those  for  pairs  5  to  9  which  are  not 
shown . 

At  this  stage  it  was  not  possible  to  take  samples 
every  revolution  and,  therefore,  the  test  time  was  much 
longer  than  expected.  This  problem  has  since  been  solved 
[10]  and  examples  are  given  below  when  measurements  from 
the  TX  are  discussed. 

In  Fig.  12  results  of  taking  paced  (synchronized) 
samples  of  a  sine-wave,  produced  by  a  signal  generator,  are 
given  for  four  different  numbers  of  repetitions.  Again,  the 


"Pacer"  was  used  as  the  trigger  device.  The  number  of  re¬ 
petitions  simply  indicates  how  often  a  sample  at  a  particu¬ 
lar  location  was  taken.  An  average  value  was  calculated 
and  is  shown  plotted  in  Fig.  12.  A  trace  of  non-uniformity 
is  seen  only  for  the  case  of  a  single  sample. 

The  PFG  was  also  used  to  obtain  comparisons  of  the 
response  of  different  kinds  of  pneumatic  and  high  response 
probes  similar  to  those  given  in  [7],  The  goal  was  to  get 
additional  information  to  add  to  the  work  reported  in  [26], 
[33]  and  [34],  for  example,  particularly  to  obtain  results 
at  higher  frequencies.  Unfortunately,  the  small  dimensions 
of  the  jet  and  wheel  did  not  allow  controlled  variations  of 
frequency  alone  to  be  made  and,  again,  the  learning  experi¬ 
ence  was  the  main  benefit  derived  from  the  experiment.  In 
Section  VI  suggestions  are  given  on  how  more  detailed  in¬ 
vestigations  could  possibly  be  carried  out  in  the  RGG  at 
the  DFVLR. 

2.  Calibration  of  New  Type  A  and  Type  B  Probes 

In  [ 3 ] ,  the  procedure  for  calibrating  the  Type  A 
and  Type  B  probes  in  a  free  jet  is  described.  The  old  probes 
in  [3]  were  the  ones  shown  in  Fig.  4.  Two  new  probes  were 
bought  which  were  smaller  in  size  in  the  sensing  area.  The 
diameter  of  the  head  was  1.63  mm  (instead  of  2.36  mm).  Addi¬ 
tionally,  the  inclination  angle  of  the  B  probe  was  35  degrees 
(instead  of  55  degrees)  as  was  used  in  [3].  The  new  probes 
were  both  derived  from  Kulite  Type  XB-062-25  transducers. 


The  probes  were  calibrated  in  the  107  mm  diameter 


free  jet  exhausting  into  the  room  and  supplied  from  the 
laboratory  compressed  air  system.  A  sketch  of  the  calibra¬ 
tion  apparatus  is  shown  in  Fig.  13.  The  probes  were  mounted 
in  a  probe  holder  which  allowed  adjustments  to  be  made  in 
yaw  and  pitch  of  the  probes  about  the  probe  tips.  The  Mach 
number  of  the  jet  was  controlled  from  approximately  .2  to  .7 
in  steps  of  .1.  The  Mach  number  was  obtained  from  the  stag¬ 
nation  pressure  measured  just  upstream  of  the  free jet  nozzle, 
and  atmospheric  static  pressure.  The  upstream  stagnation 
pressure  was  determined  to  be  equal  to  the  stagnation  pres¬ 
sure  at  the  test  probe  station  to  within  measurement  accura¬ 
cy.  All  pneumatic  pressures  were  recorded  using  a  single 
strain  gauge  transducer,  signal  conditioning  and  Scani valve, 
with  the  Steady  State  Data  Acquisition  System  shown  in  the 
upper  half  of  Fig.  7.  The  reference  sides  of  the  Kulite 
transducers  were  vented  to  atmosphere  or  supplied  with  con¬ 
trolled  air  pressure  (also  recorded  through  the  Scanivalve 
system) .  The  yaw  angle  of  the  probes  was  recorded  by  the 
data  system  from  rotational-position  potentiometers  to  an 
uncertainty  of  about  0.05°.  The  pitch  angles  of  the  probes 
were  set  manually  and  read  against  a  vernier  scale  to  an 
uncertainty  of  0.5°. 

The  calibration  was  carried  out  for  pi  ..  i  angles 
-15°  ^  ^  25°  in  steps  of  5°  and  yaw  angles  -80°ia£80°.  150 

measurements  were  taken  while  rotating  the  probe  from  a  *  8  0°  to 
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a  *  -80°.  During  this  procedure  the  probes  were  yawed  con¬ 
tinuously  but  slowly  enough  to  acquire  data  with  verified 
accuracy . 

a.  Results  for  the  Type  A  Probe.  For  an  arbitra¬ 
rily  chosen  flow  condition  the  output  of  the  probe,  p  ,  is 

A 

shown  as  a  function  of  yaw  angle,  ct,  in  Fig.  14.  In  Fig.  15 
results  for  different  Mach  numbers  are  given.  In  this  case 
data  are  plotted  only  for  values  pA  £  0  in.  H20.  In  [3]  it 
was  reported  that  for  a  constant  pitch  angle,  <f>,  the  probe 
reads  static  pressure  (indicated  by  pA  =  0)  at  approximately 
the  same  yaw  angle  position  independent  of  the  Mach  number. 
This  can  be  seen  in  Fig.  15  where  the  static  pressure  is 
registered  at  a  yaw  angle  of  a  =  ±64°  for  each  Mach  number. 
This  is  shown  in  more  detail  for  an  arbitrarily  chosen  flow 
condition  in  Fig.  16,  in  which  a  larger  ordinate  scale  has 
been  used.  At  positive  yaw  angle,  all  curves  cross  the 
abscissa  at  the  same  location.  This  is  almost  but  not 
exactly  the  situation  for  the  negative  yaw  angle. 

An  exception  to  the  preceding  statement  occurred 
for  the  lowest  Mach  number,  M  =  .2,  and  for  the  complete 
range  of  pitch  angles  investigated.  Reynolds  number  effects 
were  suspected,  but  unfortunately  this  parameter  could  not 
be  changed  independently  in  the  calibration  set-up  over  a 
sufficiently  wide  range. 

The  data  reduction  procedure  for  the  calibration 
followed,  initially,  the  method  described  in  [3].  For 
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convenience,  a  FORTRAN  computer  program  was  written  to  carry 
it  out.  The  quantities  which  must  be  evaluated  are  described 
in  Section  IV.  Here  is  described  only  how  the  coefficients 
are  determined. 

First,  the  non-dimensional  pressure  coefficient, 
c  ,  defined  in  equation  (1) ,  is  computed.  (Data  which  are 

Sr 

obviously  incorrect  are  ignored.)  The  probe  yaw  angles,  ag, 
where  static  (atmospheric)  pressure  is  indicated  are  evaluated 
by  applying  a  least  squares  linear  curve  fit  through  nine  mea¬ 
sured  values  which  are  close  to  the  expected  value  of  otg.  In 
addition  to  the  smallest  positive  cp  value  the  adjacent  4 
positive  and  4  negative  values  are  taken. 

A  small  difference  in  magnitude  was  found  in  the 

positive  and  negative  values  of  a  (up  to  2  degrees)  which 

s 

indicates  the  presence  of  a  small  error  in  setting  the  probe 

to  zero  in  the  holder  when  aligned  with  the  flow.  Therefore, 

the  average  magnitude  from  positive  and  negative  values  of 

ct  is  taken, 
s 

The  computation  of  the  maximum  c^  value,  cp0A/ 
is  done  by  using  a  least  squares  curve  fit  of  fourth  order 
through  the  measurement  points  which  are  within  the  range 
-30°  *  a  +30°.  When  the  same  method  is  applied  for  the 
complete  range  of  a  it  leads  to  deviations  from  the  data 
which  are  largest  in  the  area  of  the  maximum.  Simply  asking 
for  the  highest  actual  value  of  cp  can  result  in  a  value 
which  is  too  large  because  of  the  experimental  scatter. 
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(Only  obviously  wrong  values  are  ignored.)  Of  course  the 
maximum  should  be  expected  to  occur  at  a  =  0.  But,  again 
because  of  scatter  this  is  not  assumed.  Although  deviations 
of  the  value  of  a  corresponding  to  the  calculated  value  of 
CpgA  and  the  value  of  a  which  results  from  selecting  the  mid¬ 
point  between  positive  and  negative  values  of  a  can  be  up  to 
10  degrees,  only  negligible  differences  in  the  value  CpQ  it¬ 
self  occur  because  of  the  flatness  of  the  probe  characteristic 
near  its  area  of  maximum  output. 

As  an  additional  example  in  Fig.  17  plots  for 
Cp/cpgA  “  f(«)  are  shown  for  M  »  .6  and  25°  £  c*  1  -10°.  After 
non-dimensionalising  a  with  ctg,  54  different  curves  for 
Cp/cp0A  *  f(a/ag)  are  obtained  and  a  fourth  order  polynomial 
curve  is  fitted  by  the  least  squares  method  to  each  of  the  54 
data  sets. 

For  -1.1  -  a/a  £  +1,1  in  steps  of  .1  23  values 
of  Cp/CPQA  are  calculated  for  each  of  the  54  polynomials. 

1242  values  are  then  used  to  calculate  5  coefficients  for  a 
single  fourth  order  representative  curve.  Fig.  18.  Because 
of  the  fact  that  the  probe  behaves  somewhat  differently  at 
M  *  .2  additionally,  coefficients  are  provided  which  exclude 
the  test  at  this  Mach  number.  For  comparison  the  result  for 
the  old  A  probe  is  also  plotted. 

The  deviation  of  the  results  for  M  -  .2  is  even 
more  distinct  when  CpQA  is  plotted  versus  pitch  angle,  <f>. 

Figs.  19  and  20,  and  ag  versus  <p ,  Figs.  21  and  22.  These 
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are  additional  curves  which  are  used  in  the  reduction  of  data 
from  TX  measurements.  The  data  in  Figs.  19  and  21  are  shown 
connected  for  each  Mach  number  by  straight  lines.  For  each 
of  the  curves  a  fourth  order  curve  fit  is  computed,  and  with 
the  aid  of  the  coefficients  25  values  are  calculated  in  the 
range  of  interest.  With  the  resulting  150  numbers  two  curves 
as  described  in  the  last  paragraph  are  evaluated.  Figs.  20 
and  22. 

b.  Results  for  the  Type  B  Probe.  For  the  Type  B 
probe  only  CpgB  as  a  function  of  0  is  required  [3].  The  pro¬ 
cedure  used  for  acquiring  calibration  data  is  the  same  as  was 
described  for  the  A  probe.  Because  of  the  35°  inclination 
the  B  probe  shows  a  totally  different  characteristic  in  com¬ 
parison  to  the  A  probe.  This  can  be  seen  from  results  shown 
in  Figs.  23  and  24,  which  are  representative  examples  of  the 
acquired  data.  The  final  result  is  given  in  Fig.  25  together 
with  the  coefficients.  The  straight  lines  reqresent  the  ex¬ 
perimental  data  and  the  dotted  ones  are  the  results  after 
the  curves  are  fitted.  It  is  to  be  seen  that  a  distinct  Mach 
number  dependence  occurs.  That  was  not  the  case  for  the  old 
B  probe  and,  therefore,  only  one  curve  is  plotted  in  Fig.  25. 

From  a  first  attempt  to  calibrate  the  probe  a  se¬ 
verely  unsymmetric  characteristic  was  measured  which  was  sub¬ 
sequently  found  to  be  due  to  some  very  small  dust  particles 
in  some  holes  of  the  screen  [45].  When  the  screen  was  cleaned 
(which  had  to  be  done  extremely  carefully)  a  nearly  perfect 


behavior  resulted.  This  experience  is  mentioned  because  dust 
particles  cannot  be  avoided  completely  in  any  research  fa¬ 
cility  like  the  TX.  It  is  suggested  that,  when  data  are  sus¬ 
pect  and  instrumentation  has  been  checked,  this  occurrence  of 
fouling  should  be  taken  into  account  as  a  possible  cause  of 
incorrect  measurements.  The  question  can  be  resolved  by  a 
visual  inspection  of  the  probe  tips  under  a  microscope. 

3 .  Application  of  Reduction  Programs  to  Calibration  Data 

The  programs  which  were  written  to  evaluate  the 
velocity  vecotrs  from  probe  measurements  were  applied  first 
to  the  calibration  data.  The  idea  was  to  provide  results 
which  would  indicate  the  degree  of  accuracy  in  applying  the 
approximations.  In  [3]  a  similar  test  is  described.  There, 
both  probes  were  physically  remounted  in  the  calibration  rig 
at  corresponding  yaw  and  pitch  angles  and  verification  of 
different  flow  conditions  were  calculated  using  the  reduc¬ 
tion  program. 

In  the  case  of  the  new  probes,  remounting  and  re¬ 
testing  the  probes  was  not  necessary  since  calibration  data 
were  taken  at  entirely  similar  flow  conditions  for  the  two 
probes.  The  reduction  procedure  was  simply  applied  to  the 
data  taken  during  calibrations  and  the  results  compared  with 
the  known  conditions  set  in  the  test  apparatus. 

a.  Results  for  Analytical  Data  Reduction.  Figs . 

26  to  28  show  some  of  the  results  of  several  that  were  ob¬ 
tained  using  the  method  described  in  Section  IV.  Diagrams 
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in  Figs.  26  and  27  are  results  from  the  old  probes  and  those 
in  Fig.  28  are  from  the  new  pair  of  probes.  Using  Eq.  (6) , 
for  data  at  fixed  Mach  number  and  pitch  angle  the  values  for 
pQA  and  PqB  are  constant  and  pA(m)  and  otA(m)  are  the  probe 
pressure  and  yaw  angle.  Mach  number,  K,  and  pitch  angle, 

<j>,  can  be  evaluated  for  each  value  of  pA(m)  and  aA(m)  and 
these  are  shown  in  Figs.  26  to  28.  The  fixed  test  condition 
is  indicated  as  a  reference  value.  It  is  to  be  seen  that 
the  method  gives  unacceptably  large  errors  in  the  area  of 
aobut  -40°  <  a  <  +40°,  where  dotted  lines  are  used  to  connect 
corresponding  results.  This  is  to  be  expected.  But  even  in 
remaining  regions  the  results  are  not  accurate,  especially 
the  Mach  number  which  depends  (in  an  iterative  process)  on 
the  evaluation  of  the  static  pressure.  (The  values  for  total 
pressure  are  quite  acceptable.)  In  the  pitch  angle  there 
seems  to  be  a  better  consistency  because  deviations  from  the 
test  condiitons  are  reasonably  constant  in  the  area  where 
good  results  are  expected. 

Those  general  remarks  are  true  for  both  pairs  of 
probes.  For  the  new  probes,  much  more  calibration  data  are 
available  and,  therefore,  the  evaluated  coefficients  and 
corresponding  results  were  initially  expected  to  be  better 
than  those  for  the  old  probes .  But  of  this  there  is  no  in¬ 
dication  in  the  plots.  The  reason  which  is  suggested  at  the 
moment  is  that  the  procedure  is  sensitive  to  the  accuracy  of 
the  fourth  order  curve  fit.  The  calibration  (vs  yaw)  is 
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represented  quite  well  for  the  B  probe  and  for  the  A  probe 
away  from  the  neighborhood  of  a  **  0,  but  it  may  be  that  it  is 
not  good  enough  for  the  A  probe  near  to  a  =  0.  In  [8]  a  sine 
function  was  used  to  represent  the  calibration  data  for  an  A 
probe.  The  function  provided  a  good  fit  to  the  "flattened" 
type  of  characteristic.  It  is  possible  that  the  use  of  a 
similar  function  for  the  (smaller)  new  probes  would  improve 
the  results  obtained  with  the  analytical  method  of  reduction. 

In  Fig.  28  results  are  shown  for  reductions  car¬ 
ried  out  using  coefficients  evaluated  for  different  Mach  num¬ 
ber  ranges.  Results  for  symbols  x  and  +  are  evaluated  using 
coefficients  given  in  Figs.  18,  20,  22  and  25.  It  is  to  be 
seen  that  the  curves  always  show  similar  trends  and  are  with¬ 
out  significant  differences  in  quality  with  respect  to  accur¬ 
ately  reproducing  the  known  test  conditions.  What  is  somewhat 
surprising  is  that  even  the  data  reduced  using  the  coefficients 
evaluated  for  the  very  same  flow  Mach  number  and  pitch  angle 
condition  set  in  the  calibration  rig  (M  =  .5  and  <j>  =  0°), 
which  are  shown  plotted  with  symbol  o,  show  a  very  similar 
behavior.  Therefore,  it  is  likely  that  the  problem  lies  in 
the  particular  choice  of  analytical  function  and  the  high 
sensitivity  in  an  analytical  approach  of  this  kind. 

b.  Results  for  Numerical  Data  Reduction.  The  nu¬ 
merical  reduction  method  of  [14] — Method  II,  was  used  to  re¬ 
duce  the  calibration  data  for  the  old  A  and  B  probes.  The 
available  results  from  the  TX  were  obtained  with  this  set  of 
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probes  and  therefore  a  test  of  the  procedure  for  these  probes 
was  desirable. 

In  Figs.  29  and  30  the  results  for  pitch  angle, 

<J>,  yaw  angle,  a,  and  Mach  number,  M  are  shown  plotted  versus 
±a.  Using  the  abscissa  in  this  way  should  indicate  that  two 
(of  the  four  required) .  pressures  are  taken  from  the  left  and 
right  branches  of  the  calibration  curve  for  the  same  absolute 
value  of  a.  The  other  two  pressures  are  the  maximum  outputs 
of  each  of  the  A  and  B  probes.  It  should  be  mentioned  that 
every  other  combination  of  pressures  (which,  in  principle,  it 
is  possible  to  use)  gave  poorer  results  in  comparison  to  the 
aforementioned  choice.  For  example,  it  is  possible  to  use 
three  values  from  the  B  probe  and  only  one  from  the  A  probe, 
and  even  four  values  from  only  one  probe  will  lead  to  a  suc¬ 
cessful  solution  of  the  iteration.  It  is  stated  however  in 
[14],  and  now  confirmed  by  experience,  that  using  both  maxi¬ 
mum  values  of  the  probes  and  two  others  which  correspond  to 
symmetric  yaw  angle  positions  will  give  the  best  results. 
Indeed,  the  yaw  position  which  is  recommended  in  [14]  is 
±25°,  but  for  real  values  of  calibration  data  yaw  angles  in 
the  range  of  about  50°  to  65°  must  be  used.  This  is  quite 
clearly  indicated  in  the  diagrams  of  Figs.  29  and  30.  It  is 
to  be  seen  that  in  this  range  of  a  the  required  flow  quanti¬ 
ties  show  an  asymptotic  behavior  (for  simplification  the  re¬ 
sults  are  connected  by  straight  lines)  and  the  results  are  in 
quite  reasonable  agreement  with  the  test  conditions.  In 


areas  near  maximum  probe  outputs  (±a  <  25°)  incomprehensible 
results  occur  as  was  the  case  of  the  analytical  method 
(-40°  <  a  <  +40°) . 

The  scatter  of  the  results  is  less  in  comparison 
to  that  found  using  the  analytical  method  (Figs.  26  to  28, 
where  the  same  scales  are  used) .  Therefore,  for  the  moment 
the  numerical  method  is  recommended  as  the  more  successful 
one. 

If  ±ot  values  are  used  which  are  lower  than  those 
from  the  recommended  range,  yaw  angle  deviations  from  the 
test  conditions  increase  very  rapidly.  The  behavior  of  the 
yaw  angle  calculated  using  data  at  large  yaw  angles  is  not 
as  obviously  asymptotic  as  are  the  Mach  number  and  pitch 
angle.  This  can  be  interpreted  as  an  indication  of  the  sen¬ 
sitivity  of  the  reduction  procedure.  In  view  of  the  magni¬ 
tude  of  the  possible  error  in  a  it  might  be  better  to  derive 
the  yaw  angle  as  described  for  the  analytical  method  in  an 
earlier  section  from  the  position  of  the  maximum  of  the  probe 
output.  It  is  suggested  that  these  different  results  for 
yaw  angles  be  compared.  This  topic  will  be  discussed  in 
connection  with  the  results  from  the  TX. 

B.  RESULTS  FROM  THE  TX 

Results  from  high  response  pressure  transducers  are  ob¬ 
tained  with  the  A  and  B  probe  system  in  the  measurement 
plane  between  rotor  and  stator,  and  with  twelve  different 
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labeled  K6  to  K13  a  number  (see  Fig.  5) .  The  number  indi¬ 


cates  the  axial  location,  where  9  is  at  the  rotor  entrance 


plane  and  smaller  numbers  indicate  planes  upstream  from 
there. 

An  example  of  the  data  obtained  from  the  Type  A  and  B 
probes  using  the  acquisition  procedure  described  in  Section 
III  are  given  in  Fig.  31.  Some  general  remarks  can  be  made. 
Outputs  for  both  probes  are  shown  plotted  for  one  blade  pair 
for  two  different  numbers  of  repetitions  (compare  Fig.  12) . 
[Note  that  all  results  for  the  A  and  B  probes  presented  in 
this  report  are  taken  at  the  mid-span  position  (50%  of  blade 
height) . ]  The  number  of  counts  in  the  circumferential  direc 
tion  is  indicated  by  n;  one  passage  is  divided  into  128 
counts.  The  noise  reduction  in  fluctuations  which  accom¬ 
panies  increasing  the  number  of  repetitions  is  clearly  seen. 
Of  course,  this  would  improve  if  even  higher  numbers  of 
repetitions  were  used.  With  the  software  improvements  de¬ 
scribed  in  [10],  and  with  the  number  of  repetitions  equal 
to  ten,  the  data  acquisition  time  has  been  reduced  to  about 
10  secs  from  the  45  seconds  which  were  needed  with  the  old 
software  in  [6].  Additionally,  it  should  be  mentioned  that 
the  averaged  value  indicated  in  the  diagrams  by  a  broken 
line  is  nearly  independent  of  the  number  of  repetitions 
used.  The  wakes  of  the  blades  are  clearly  indicated  and  it 
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is  to  be  seen  that  they  occur  almost  exactly  128  counts 
apart,  which  is  another  indication  of  the  precision  of  the 
high  speed  data  acquisition. 

In  Fig.  32  pictures  taken  from  an  oscilloscope  are  pre¬ 
sented  to  demonstrate  transducer  outputs  for  the  A  and  B 
probe  and  for  the  wall  statics,  in  this  particular  case  K9.5 
and  K10.  It  is  clear  from  the  l/30th  sec.  exposure  that 
much  more  noise  occurs  for  the  probes  inserted  into  the  flow 
than  for  the  measurements  made  at  the  wall  even  if  the  dif¬ 
ferent  scales  are  taken  into  account.  However,  results  in 
Fig.  31  also  indicate  that  increasing  numbers  of  repetitions 
will  result  in  well-defined  average  signals.  This  would  not 
at  first  be  expected  from  watching  the  probe  outputs  on  the 
oscilloscope. 

Results  for  A  and  B  probes  can  now  be  discussed  in  more 
detail. 

1.  Results  for  A  and  B  Probes 

In  Figs.  33  and  34,  results  for  different  numbers  of 
repetitions  are  shown  (in  this  case  for  the  B  probe  only) . 

In  these  diagrams,  following  the  DPDS  procedure  described  in 
Section  IV,  signals  are  shown  from  the  B  probe  when  set  at 
specific  yaw  angles,  a  (for  definition  see  Fig.  3) .  It  is 
noted  that  the  averaged  values  for  corresponding  curves  with 
different  numbers  of  repetitions  would  differ  slightly,  be¬ 
cause  data  were  taken  on  different  days.  It  is  to  be  seen 
that  with  increasing  values  for  a  the  basic  behavior  of  the 


signals  within  the  wakes  changes.  For  example,  in  the  case 
of  a  =»  -43°  within  the  wakes  extrema  are  indicated  as  minima 
and  in  the  case  of  0(3+103°,  maxima  occur.  The  difference 
between  peak  to  peak  output  values  decreases  for  a-values 
which  are  between  these  extremes. 

With  the  aid  of  Fig.  35  a  qualitative  explanation 
of  this  behavior  can  be  given.  Here  velocity  vectors,  probe 
positions  and  curves  of  "probe  output,  p,  versus  angle,  a," 
are  sketched.  (It  is  of  no  importance  for  the  following 
qualitative  description  but  it  should  be  mentioned  to  avoid 
misunderstanding  that  these  curves  are  actually  different 
for  different  positions  of  the  probe  across  the  blade  pas¬ 
sage.)  If  the  rotor  blade  wake  is  a  simple  wake,  the  rela¬ 
tive  velocity  in  the  wake,  w2D,  is  smaller  than  that  outside 
the  wake,  w2,  but  the  relative  flow  angle  is  almost  constant 
This  can  be  seen  to  increase  the  absolute  flow  angles  seen 
by  the  probe  (a2  *  outside  the  wake;  a2D  35  inside  the  wake)  . 
This  results  in  probe  outputs  pc2  and  pc2D  for  a  “  -43°  as 
sketched  on  the  right  hand  side  of  Fig.  35.  Here,  PC2  > PC2D 
but  both  have  negative  values  (that  means  below  atmospheric 
pressure),  compare  the  result  for  a  =*  -43°  in  Fig.  34.  If 
the  probe  angle,  a,  is  increased  both  values  pc2  and  pc2D 
are  "traveling"  on  curves  p  ■  f(a)  in  the  positive  a  direc¬ 
tion.  For  example,  in  the  case  of  a  *  +52°  both  values  are 
quite  the  same,  and  therefore  distinct  peak-to-peak  differ¬ 
ences  do  not  occur.  For  a  =  103°  Pc2D  >  Pc2  which  explains 
the  "inverted"  signal  behavior. 
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It  is  recommended  in  [3]  that  measurements  for  seven 
different  yaw  angle  positions  should  be  obtained  to  enable 
an  adequate  reduction  of  data  from  the  DPDS  technique.  The 
first  data  were  obtained  (examples  are  shown  in  Figs.  33  and 
34)  following  this  suggestion.  However,  it  was  subsequently 
found  that  more  angle  settings  were  required  to  improve  the 
data  reduction  procedure.  Both  acquisition  program  [16] 
and  reduction  program  allow  measurements  at  eleven  different 
probe  settings,  and  this  feature  was  used  in  subsequent 
tests. 

Two  sets  of  data  from  the  TX  using  the  old  A  and  B 
probes  are  shown  in  Figs.  36  to  39.  They  can  serve  to  il¬ 
lustrate  some  general  statements  with  respect  to  the  DPDS 
technique.  (Following  the  failure  of  the  A  probe  as  a  re¬ 
sult  of  handling,  further  measurements  were  not  possible 
with  this  system.) 

Each  section  of  Figs.  36  to  39  shows  results  for 
both  A  and  B  probes  for  different  blade-to-blade  positions 
and  each  represents  a  discrete  value  of  the  probe  angle 
setting.  For  the  first  data  set  in  Fig.  36  and  Fig.  37  the 
number  of  repetitions  was  10  and  for  the  second  set  in  Fig. 
38  and  Fig.  39,  40  repetitions  were  used. 

Comparing  the  corresponding  signals  from  the  two 
data  sets,  it  can  be  seen  that  for  the  first  set  an  apparent 
"phase"  shift  occurs  between  the  A  and  B  probe  results  which 
is  not  the  case  for  the  second  data  set.  In  fact  a  change 
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was  made  in  the  data  acquisition  software  following  the  first 
data  set.  It  was  found  that  one  probe  was  placed  at  a  dif¬ 
ferent  circumferential  location  (135°  instead  of  90°,  Fig. 

5)  after  a  bearing  failure  of  the  drive  turbine  was  repaired 
and  the  instrumentation  was  replaced.  Therefore,  the  number 
of  counts  separating  the  two  probes  which  was  programmed  in 
the  acquisition  software  (to  provide  measurements  with  the 
two  probes  from  the  same  blade  pair)  did  not  match  the  hard¬ 
ware  any  more.  After  correcting  the  number  of  counts  in  the 
program,  the  signals  became  aligned  correctly,  as  seen  in 
Figs.  38  and  39. 

There  is,  indeed,  no  obvious  check  which  would  en¬ 
sure  that  measurements  are  obtained  from  the  same  blade 
pairs,  but  the  test  with  the  PFG  described  in  Section  V.A.l 
indicated  clearly  that  the  trigger  system  worked  correctly, 
when  correctly  programmed. 

The  phase-shift  of  the  first  set  of  data  could  be 
corrected  partially  by  providing  an  offset  of  a  fixed  number 
of  counts  in  the  reduction  program.  28  counts  were  needed, 
but  the  disadvantage  remains  that  signals  were  taken  behind 
different  blade  pairs.  Furthermore,  a  complete  reduction 
over  two  passages  was  not  possible,  unless  the  results  which 
were  offset  to  negative  counts  were  patched  to  the  data  at 
the  other  end  to  complete  the  two  passages  (a  procedure 
which  may  be  doubtful) . 
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The  corrected  software  in  [10]  needed  test  times  of 
about  25  minutes  for  one  complete  set  of  data  from  2  probes 
at  a  fixed  radial  location,  taken  at  eleven  probe  settings 
with  40  repetitions.  A  third  of  this  time  was  required  to 
manually  adjust  the  probe  yaw  angles.  Although  the  ratio 
of  "test  time  to  quality  of  signal"  is  acceptable,  it  may  be 
worthwhile  to  provide  remotely  actuated  probes  in  the  future 
firstly  to  decrease  the  test  time  and  secondly  for  the  con¬ 
venience  and  safety  of  the  operator. 

Applying  the  data  reduction  (with  on-line  calibra¬ 
tion)  to  the  signals  obtained  from  the  A  and  B  probes  re¬ 
sults  in  2  x  256  different  curves  of  p  =  f(a)  as  shown  for 
the  first  data  set  in  Figs.  40  to  45.  The  curves  shown  are 
the  result  of  plotting  probe  pressure,  p,  for  different  cir¬ 
cumferential  locations,  n  (counts)  for  the  11  yaw  angles,  a. 
Measured  data  points  are  connected  with  straight  lines  for 
clarity.  In  Figs.  40  to  45  the  curves  are  drawn  to  a  scale 
which  is  shown  once  at  the  bottom  for  the  lowest  curve  and 
the  other  curves  each  offset  in  sequence  by  about  3.4  inches 
of  water.  In  Figs.  42  and  45  the  results  in  the  area  of  the 
wake  of  the  first  blade  is  shown  in  more  detail,  with  each 
curve  offset  by  almost  seven  inches  of  water. 

A  second  set  of  data  were  obtained  at  a  highly  throt' 
tied  flow  condition  where  a  higher  degree  of  flow  separation 
can  be  expected  to  occur  in  the  rotor  blading.  The  same  yaw 
angle  settings  were  used.  Plots  of  p  *  f(a)  for  these  data 


indicate  that  especially  within  the  wakes  Cparticularly  for 
the  A  probe)  yaw  angles  greater  than  103°  were  needed  in 
order  to  fuuly  describe  the  right  side  of  the  characteristic 
curve.  An  example  is  shown  in  Fig.  46.  This  is  a  diffi¬ 
culty  because  the  reduction  procedure  requires  adequate  data 
on  both  left  and  right  branches . 

In  wake  flow  areas  the  probe  signal  seems  to  be 
relatively  more  unsteady  when  the  probes  are  set  to  flow 
angles  which  are  close  to  the  mean  flow  direction.  This  can 
be  seen  in  the  A  probe  signals  shown  in  Figs.  46a  and  b. 

For  a  probe  angle  which  is  close  to  the  aligned  flow  direc¬ 
tion  (a  =  52°,  Fig.  46a)  it  is  seen  that  within  the  wake 
the  noise  is  relatively  more  pronounced  than  it  is  outside 
the  wake.  For  the  same  numbers  of  repetitions,  at  large 
values  of  probe  yaw  angles  (a  »  103°,  Fig.  46b) ,  the  noise 
inside  the  wake  flow  appears  to  be  very  similar  to  that  out¬ 
side  the  wake  flow  and  quite  similar  to  that  outside  the 
wake  flow  at  the  probe  angle  of  52°. 

In  view  of  this  there  are  two  possible  ways  to  in¬ 
crease  the  quality  of  the  signals  without  extending  the  data 
acquisition  time.  First,  the  number  of  repetitions  can  be 
increased  for  probe  angles  which  are  close  to  flow-aligned 
conditions  and  can  be  decreased  for  the  rest  of  the  angles; 
or  second,  for  all  probe  angle  settings  the  number  of  repe¬ 
titions  can  be  increased  when  data  are  acquired  in  wake  flow 
areas  and  can  be  decreased  outside  those  areas.  The 


identification  of  wake  flow  areas  can  be  made  from  a  display 
of  the  transducer  output  on  the  screen  of  the  oscilloscope. 

At  first  glance,  a  comparison  of  the  curves  shown 
in  Figs.  40  to  46  with  those  from  the  calibration,  Fig  17 
for  example,  will  lead  to  the  conclusion  that  the  test  data 
result  in  pressure  characteristic  curves  similar  to  those 
obtained  in  calibration  tests.  And  it  may  be  considered  sur¬ 
prising  that  signals  like  those  shown  for  example  in  Figs. 

36  and  37  result  in  quite  well-shaped  curves,  especially 
outside  the  wakes;  moreso,  in  view  of  the  unsteadiness  evi¬ 
dent  in  the  oscilloscope  traces  (Fig.  32).  Nevertheless, 
going  further  into  the  details  it  turns  out  that  the  initial 
statement  does  not  hold  true  everywhere,  especially  on  the 
ends  of  the  left  and  right  branches  and  particularly  within 
the  wake  flow  areas.  It  is  known  from  applying  the  reduction 
methods  to  the  calibration  data  that  pressures  must  be  used 
which  are  obtained  at  probe  yaw  angles  which  are  approxi¬ 
mately  ±60°  away  from  the  angle  for  maximum  probe  output. 
Techniques  are  discussed  below  which  may  help  to  overcome 
this  very  important  problem. 

While  being  aware  of  the  aforementioned  problem  in 
practice  a  data  reduction  procedure  is  described  next  which 
is  based  on  the  assumption  of  having  a  perfect  pressure 
characteristic  curve.  This  is  done  to  gain  a  better  under¬ 
standing  and  illustrate  the  initial  problems  of  the  DPDS 
technique. 


From  the  first  data  set,  the  results  at  n  *  31  Cor 
n  *  3,  if  the  added  offset  of  28  counts  is  not  included)  are 
shown  plotted  in  more  detail  in  Fig.  47.  It  is  noted  that 
it  was  carefully  verified  that  applying  the  offset  did  not 
change  the  shape  of  the  probe  characteristic. 

Using  a  fourth  order  curve-fit  to  the  data  points 
shown  plotted  with  the  symbol  •  (through  which  a  curve  has 
been  hand-faired)  and  calculating  corresponding  points  using 
the  fourth  order  polynomial  results  in  pressure  values  indi¬ 
cated  by  the  x  symbols.  It  is  seen  that  in  this  case  the 
polynomial  curve  fit  was  a  good  representation  of  experi¬ 
mental  results.  Next,  a  line  is  drawn  parallel  to  the  ab¬ 
scissa  and  moved  with  respect  to  the  abscissa  to  find 
iteratively  the  pressure  p  ,  where  the  length  of  the  line 
between  the  sides  of  the  curve  is  equal  to  128°.  This  pro¬ 
cedure  is  used  because  calibration  data  showed  that  static 
pressure  is  sensed  by  the  A  probe  at  about  64°  (or  ±64°) 
independently  of  Mach  number  and  pitch  angle  (Figs.  16  and 
22).  In  Fig.  47  values  are  indicated  for  p„  and  aT ,  a_  and 
the  corresponding  midpoint,  <*A,  which  is  regarded  as  the 
flow  yaw  angle. 

In  addition,  it  is  assumed  that  the  A  probe  indi¬ 
cates  total  pressure,  pfc,  in  its  maximum  output  value,  which 
is  very  reasonable  as  long  as  pitch  angles,  $,  are  roughly 
in  the  range  j $ |  <  12°.  In  the  ideal  case,  the  maximum 
pressure,  p  ,  which  is  evaluated  from  the  fitted  curve, 


should  be  located  at  aA,  but  a  lack  of  symmetry  in  the  charac 
teristics  is  known  to  be  present  from  the  calibration  tests. 
The  difference  between  otA  and  the  angle  for  maximum  output  is 
about  7  degrees  in  the  example  shown.  Because  of  the  flat 
behavior  in  the  neighborhood  of  the  curve's  maximum,  p(aA)  is 
used  as  total  pressure,  pt,  instead  of  Pmax*  The  ratio  Ps/Pt 
then  allows  the  calculation  of  the  local  Mach  number,  M. 

The  assumption  that  the  A  probe  is  not  indicating 

static  pressure,  pg,  at  aR  -  aL  *  128°  but  ±1°  from  this 

value  will  lead  to  deviations  in  Mach  number  of  about  ±2% 

from  the  previously  calculated  value.  The  general  trends  of 

the  calculated  Mach  number  distribution  for  a  set  of  data 

would  not  change  but  the  average  value  would  change.  The 

same  situation  will  occur  if  the  yaw  angle  is  evaluated  at 

p  and  pm,„  is  taken  to  be  the  stagnation  pressure, 
max  ^max 

For  both  Mach  number,  M,  and  yaw  angle,  a  results  of 
the  aforementioned  procedure  for  all  n  -  256  counts  are  shown 
in  Fig.  48  for  the  first  set  of  data.  The  results  follow  a 
credible  trend,  but  the  precision  will  be  questioned  in  the 
following  discussion. 

In  principle,  the  same  procedure  can  be  applied  also 
to  the  second  data  set  (at  reduced  flow,  throttled  compressor 
conditions) .  Some  mathematical  tricks  are  needed  to  overcome 
the  problem  of  missing  data  in  wake  areas.  This  has  been 
done?  however  the  results  are  not  presented  because  the  ac¬ 
curacy  is  seriously  questioned.  In  comparison  to  the  first 


results  a  smoother  Mach  number  distribution,  higher  yaw 
angles  and  larger  wake  widths  are  calculated,  all  of  which 
could  be  interpreted  as  the  result  of  a  more  severe  flow 
separation  due  to  the  highly  throttled  test  conditions. 

The  yaw  angle  distribution  can  also  be  evaluated 
using  results  from  the  B-probe.  The  characteristic  of  this 
probe  is  slightly  unsymmetric  as  is  the  case  for  the  A  probe 
and  therefore  results  for  ctg  will  again  depend  on  the  pres¬ 
sure  level  at  which  the  calculation  is  made.  Because  of  its 
extreme  inclination  to  the  flow  even  at  zero  pitch  angle 
static  pressure  is  not  well  defined  at  a  particular  ot-value 
and  must  be  based  on  the  results  from  the  A  probe.  It  was 
found  that  a  difference  of  about  10  degrees  (dependent  on 
which  level  of  pressure  was  chosen  at  which  to  evaluate  it) 
occurred  in  the  values  of  a  calculated  using  the  A  or  the  B 
probes.  This  was  ^true  for  both  sets  of  data.  If  all  uncer¬ 
tainties  in  the  various  methods  of  evaluating  a  were  con¬ 
sidered,  the  difference  (which  was  consistently  there)  of 
10  degrees  could  not  be  explained.  However,  a  misalignment 
between  one  or  both  probes  and  the  angle  scale  on  the  probe¬ 
holders  is  a  probable  explanation  since  no  verification  of 
the  angle  zero  had  been  made  after  several  months  of  opera¬ 
tion.  Unfortunately  a  failure  occurred  in  the  A  probe  as 
a  zero  check  was  being  prepared.  Since  the  old  A  and  B 
probe  set  would  not  be  used  further,  reasons  for  the  de¬ 
viations  in  yaw  angles  were  not  resolved. 


Results  calculated  here  £or  the  flow  quantities  are 
not  considered  to  be  satisfactory  using  either  the  analytical 
or  the  numerical  data  reduction  methods,  particularly  for 
the  pitch  angle,  <J>.  It  was  first  thought  that  the  discrepan¬ 
cies  in  the  data-sets  described  earlier  (the  presence  of  an 
offset  between  A  and  B  probe  data,  and  incomplete  data  for 
the  right  side  of  the  pressure  characteristic)  might  be  re¬ 
sponsible.  Therefore,  parts  of  another  data-set  obtained  by 
F.  Neuhoff  (and  reduced  using  a  third  data  reduction  proce¬ 
dure  [45])  was  used  to  check  the  general  behavior  of  the 
characteristic  pressure  curves  derived  from  the  probe  out¬ 
puts.  This  particular  set  of  data  was  entirely  free  of  the 
disadvantages  mentioned  earlier:  data  were  taken  behind  the 
same  blade  pair,  40  repetitions  were  used  and  the  yaw  angles 
covered  the  required  range,  even  in  wake  areas.  It  was  found 
that  the  curves  of  p  =  f(a)  for  the  third  data  set  were  quite 
similar  with  respect  to  their  basic  shape  to  those  shown  in 
Fig.  40.  Consequently,  applying  either  reduction  method  to 
this  set  of  data  will  lead  to  similarly  doubtful  results. 

In  can  be  stated  therefore  that  the  basic  problem  in 
the  data  reduction  so  far  is  that  the  shape  of  both  left  and 
right  branches  of  the  characteristic  pressure  curves  ob¬ 
tained  from  test  data  do  not  fully  agree  with  the  behavior 
obtained  during  the  calibration  of  the  probes.  Since  only 
the  maximum  pressure  value  is  taken  from  the  B  probe  the 
problem  lies  only  with  the  A  probe. 
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Because  of  varying  pressure  levels  accompanying  the 
varying  Mach  numbers  and  pitch  angles  which  are  to  be  deter¬ 
mined,  a  comparison  of  the  calibration  curves  with  the  appli 
cation  curves  on  the  basis  of  raw  pressure  is  not  useful. 

For  a  more  reliable  check  the  non-dimensional  pressure  coef¬ 
ficient,  cpA#  given  by  Eq.  (1) ,  is  examined  as  a  function  of 
a.  For  the  calibration  data,  the  plots  of  cpA/cp0A  *  f(a) 
are  shown  in  Fig.  17.  Although  the  curves  are  quite  inde¬ 
pendent  of  the  different  flow  conditions,  and  the  value  of 
CpQA  depends  only  on  pitch  angle,  <t>,  and  is  very  close  to 
unity  for  0°  <  <t>  <  15°,  it  is  not  possible  to  generate  a 
direct  comparison  with  the  test  data  unless  some  assumption 
is  made  concerning  the  value  of  p_  or  of  Mach  number. 

The  total  pressure,  pt,  can  be  evaluated  reliably 

using  the  procedure  described  in  Fig.  47  but  the  static 

pressure,  p_,  in  Eq.  (1)  is  not  known.  Therefore  the  ex- 
s 

pected  Mach  number  range  was  used  as  a  parameter  and  cor¬ 
responding  p3  values  were  used  in  the  computation  of  CpA. 
This  procedure  was  applied  to  the  data  from  about  70  counts 
of  the  third  data  set  [45],  which  were  taken  at  9  different 
yaw  angles,  a. 

In  Fig.  49  one  example  is  given  for  data  obtained 
outside  the  wake  flow.  This  flow  condition  could  corres¬ 
pond,  for  example,  to  n  »  1  in  Fig.  40.  The  calibration 
curve  CpA  =»  f(a),  plotted  on  transparent  paper,  can  be  used 
to  search  "by  hand"  for  the  particular  curve  to  which  it 


best  fits.  For  the  example  shown  this  is  the  curve  for 
M  *  .42.  The  flow  angle,  a,  can  then  be  "corrected"  in  areas 
where  the  two  curves  do  not  fit.  The  resulting  values  are 
given  in  the  diagram.  It  can  be  seen  in  Fig.  49  that  the 
two  smallest  and  the  largest  values  of  aA  have  to  be  "cor¬ 
rected"  and  that  the  "correction"  for  the  smallest  value  is 
the  largest.  The  curve  containing  the  new  values  of  aA  is 
shown,  and  it  is  seen  that  it  agrees  very  well  with  the  cali¬ 
bration  curve. 

In  Fig.  50,  as  a  second  example,  results  are  shown 
for  data  obtained  within  the  wake  flow  area  (for  example  at 
n  *  43  in  Fig.  40) .  It  is  clear  that  there  is  agreement  be¬ 
tween  the  test  data  and  the  calibration  curves  at  a  higher 
Mach  number,  M  »  .48,  and  to  "correct"  aA  on  the  left  branch 
requires  corrections  which  are  about  twice  as  large  as  were 
required  in  the  first  example.  For  the  right  branch  approxi¬ 
mately  the  same  magnitude  of  correction  is  required. 

This  behavior  can  be  explained  with  reference,  again, 
to  Fig.  35.  The  value  of  the  absolute  velocity  in  the  wake 
may  increase  due  to  the  fact  that  the  relative  velocity,  w2D, 
is  reduced.  (Note  that  if  w2D -*■(),  c2d^u’^  Moreover,  the 
sketches  indicate  that  for  the  left  branch  (for  example, 
a  ■  -43°)  the  flow  direction  to  the  sensor  of  the  probe  for 
c2D  will  be  at  a  more  extreme  yaw  angle  than  for  c2  which 
may  explain  the  necessity  for  the  large  corrections.  For 
the  right  branch  the  flow  yaw  angle  is  more  extreme  for  c2 


than  for  c2ly  but  possibly  not  as  extreme  as  for  c2D  t 

a  *  -43°,  which  can  explain  the  use  of  about  the  same 
correction. 

It  should  be  mentioned  again  that  the  general  behav¬ 
ior  demonstrated  for  these  examples  is  typical  of  the  passage 
positions  from  which  they  were  taken.  A  similar  behavior 
was  found  for  both  data  sets  shown  in  Figs.  36  to  39.  The 
implication  is  that  the  data  reduction  could  be  carried  out 
with  data  from  fewer  repetitions  (10  and  40,  respectively, 
have  been  used)  if  an  explanation  (and  a  proper  correction) 
for  the  yaw  angle  problem  can  be  found.  Also  an  "incorrect" 
behavior  of  probe  pressure  with  yaw  angle  was  observed  for 
the  B  probe  measurements.  However,  a  "correction"  procedure 
for  the  B  probe  requires  an  additional  step,  since  cpB  de¬ 
pends  on  <)>,  as  can  be  seen  from  Fig.  24. 

There  are  many  possible  reasons  for  the  observed  be¬ 
havior.  What  can  be  said  at  the  moment  however  is  specula¬ 
tive  because  detailed  investigations  have  not  been  carried 
out.  For  example,  the  exact  position  of  the  probe  tip  in 
the  flow  during  rotation  must  be  questioned.  It  is  diffi¬ 
cult  to  ensure  precise  mechanical  rotation  about  the  center 
of  the  sensor  tip  when  the  probe  is  yawed  in  its  eccentric 
mount.  If  the  probe  tip  is  slightly  long,  there  is  an  error 
in  position  which  increases  with  increasing  probe  yaw  angle 
(in  both  positive  and  negative  directions) .  Different  probe 
mounts  are  used  for  calibration  compared  to  compressor 
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applications  and  while  the  probes  are  rigid  in  construction, 
small  differences  in  bending  could  occur  in  the  two  situa¬ 
tions.  It  is  recalled  that  there  was  no  opportunity  to  re¬ 
calibrate  the  first  set  of  probes.  Although  the  dimensions 
of  the  old  probes  were  as  small  as  possible  and  yet  accom¬ 
modate  the  transducer,  they  were  possibly  too  large  for  their 
intended  application.  In  Fig.  48  the  relative  size  (tip 
diameter)  of  the  old  probes  is  indicated  at  an  arbitrarily 
chosen  circumferential  location.  The  further  application  of 
the  smaller  probes  will  give  some  information  concerning 
this  problem. 

Errors  in  tip  position  on  rotation  would  result  in 
an  apparent  phase  shift  of  the  signals  obtained  at  different 
yaw  angles.  A  wake  decay  due  to  probe  movement  to  different 
axial  distances  is  demonstrated,  for  example  in  [31].  An 
apparent  decay  is  to  be  seen,  for  example,  in  Fig.  38  where 
the  positions  for  minimum  output  are  slightly  displaced  for 
the  A  and  B  probe  results  for  yaw  angles  of  a  ^  -15°.  More¬ 
over,  a  comparison  of  the  position  of  minimum  output  for 
a  =  -42°  and  a  =  -15°  for  the  A  probe  shows  a  shift  of  about 
10  counts.  However,  it  must  be  remembered  that  in  the  wake 
all  parameters  of  the  flow  are  changing  and  the  apparent 
shift  in  minima  may  be  qualitatively  correct. 

It  is  not  possible  at  the  moment  to  give  a  full  ex¬ 
planation  for  the  behavior  of  the  probes  in  the  compressor 
flow  field.  If  the  potential  explanations  speculated  above 
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are  eliminated,  then  the  general  question  of  how  high  re¬ 
sponse  pressure  probes  behave  for  different  yaw  positions  in 
unsteady  wake  flows  must  be  examined.  So  far,  other  investi 
gations  on  this  topic  are  not  known  to  the  author. 

It  should  be  mentioned  that  this  question  is  also 
particularly  relevant  when  multiple  sensor  probes  are  used. 
In  such  probes,  the  transducers  are  placed  in  various  sur¬ 
faces  away  from  the  probe's  face.  The  different  transducers 
are  differently  oriented  with  respect  to  the  flow  direction 
at  all  times,  so  that  any  difference  in  response  which  de¬ 
pends  on  orientation  will  be  difficult  to  account  for.  From 

Figs.  49  and  50  it  is  evident  here  that  at  about  |a|  >  40° 

« 

departure  from  the  calibration  behavior  begins.  Since  mul¬ 
tiple  sensor  probes  are  larger  in  size  in  comparison  to  the 
probes  described  here  it  is  suggested  that  their  unsteady 
response  will  depart  from  their  steady  behavior. 

The  numerical  data  reduction  has  been  applied  to  the 
results  modified  as  in  Fig.  49,  since  it  was  shown  in  Sec¬ 
tion  V  that  this  method  was  more  reliable  in  comparison  to 
the  analytical  one  (Figs.  26  to  30) .  The  results  are  indi¬ 
cated  with  if  x  symbol  in  Fig.  51.  Comparing  the  curves  with 
those  evaluated  for  calibration  data  (for  example  Fig.  29) , 
it  can  be  seen  that  approximately  the  same  tendency  occurs 
for  the  flow  quantities  <j>,  a  and  M.  The  results  given  for 
yaw  angle,  a,  are  here  to  be  interpreted  as  deviations  from 
the  angle,  aA,  evaluated  with  the  procedure  demonstrated  in 
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Fig.  47.  As  before,  the  best  results  are  obtained  in  the 
range  of  about  5Q°  i  ]aj  i  65°.  At  this  point,  the  data  re¬ 
duction  has  been  carried  out  as  usual  except  that  the  values 
of  a  at  the  right  and  left  branches  were  corrected.  Thus 
the  values  of  Cp  can  be  calculated  from  pfc  and  pg  for  every 
data  point  and  those  results  can  be  compared  with  the  cali¬ 
bration  curve  CpA  =  f(a).  When  this  is  done,  it  is  found, 
as  might  be  expected,  that  deviations  occur  for  about 
| a ]  <  50°.  Since  the  calibration  table  covers  only  the 
range  of  -70°  ^  a  £  +70°  deviations  from  asymptotic  values 
in  areas  |a|  >  65°  might  also  be  expected. 

When  the  method  described  in  [45]  is  applied  to  these 
data  nearly  precisely  the  same  results  are  found. 

The  value  obtained  for  pitch  angle,  <p,  was  about 
-2.5°.  In  [17]  it  is  reported  that  the  time  averaged  value 
calculated  from  measurements  with  the  combination  probe  was 
about  <}>  =  +4.5°  and,  therefore,  results  for  <J>  within  the 
wake  must  increase  considerably  to  establish  an  average 
value  of  about  the  measured  magnitude.  The  reduction  of  the 
example  given  in  Fig.  50  showed  increases  in  the  pitch  angle 
in  the  wake,  but  of  insufficient  magnitude  to  produce  an 
average  of  4>  *  +4.5°. 

This,  in  turn,  raises  the  question  of  how  reliable 
are  the  pneumatic  measurements,  in  particular  for  measuring 
pitch  angles.  This  topic  was  discussed  to  some  degree  in 
Section  IV.  A  distinct  problem  in  the  DPDS  technique  is  the 


way  the  OLC  must  be  applied  to  the  B  probe.  In  order  to  find 
out  what  would  happen  if  the  reference  pressures  of  the  B 
probe,  Pref  3  (or  equivalently  the  intercept)  were  wrong,  the 
results  were  re-evaluated  for  differences  of  +2  and  -2  inches 
H2O  respectively.  From  Fig.  51  it  can  be  seen  that  this  is  a 
relatively  minor  influence  on  both  M  and  a  but  it  causes  a 
change  in  $  of  about  ±1.5°.  In  addition,  data  were  reduced 
for  pA  values  which  were  increased  by  +2  inches  H20  when  the 
original  value  of  pB  was  used.  Results  were  equivalent  to 
the  case  of  -2  inches  H20  on  the  B  probe  and  the  original 
value  of  pA  (+  symbols  are  used  for  both  cases) . 

The  minor  influence  of  these  changes  on  a  and  M  can 
be  explained  with  reference  to  Fig.  47  although  a  different 
reduction  process  is  used.  A  different  reference  pressure 
on  the  probes  will  change  only  the  level  of  the  curve  and 
the  evaluated  yaw  angle,  aA,  would  not  change.  With  respect 
to  Mach  number,  an  impossibly  large  change  in  pressure  level 
of  say  +20  inches  H20  (pfc  »  462.4  in  H20  and  pg  *  427.4  in 
H20)  would  change  M  by  about  2%  only.  Such  a  change  of 
total  pressure  would  however  result  in  extremely  large 
changes  in  loss  coefficients. 

2.  Results  for  Wall  Static  Pressure  Transducers 

Very  briefly  some  results  are  reported  which  were 
obtained  with  12  Kulite  pressure  transducers  mounted  flush 
at  the  outer  casing  wall.  Their  peripheral  and  axial  loca¬ 
tions  are  indicated  in  Figs.  1  and  5.  Similar  measurements 
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from  the  TX  are  reported  in  [5],  [12]  and  [21]  where  8  trans¬ 
ducers  were  used.  The  programs  to  reduce  data  to  contour 
maps  used  in  [12]  and  [21],  respectively,  were  written  in 
BASIC  language.  A  FORTRAN  program  to  include  the  12  differ¬ 
ent  signals  now  available  is  in  preparation.  Such  a  program 
for  about  40  transducers  is  reported  in  [32]. 

In  52  a  series  of  signals  is  shown.  The  data 

acquisition  time  was  about  10  minutes  including  measurements 
of  all  pneumatic  pressures  and  two  OLCs.  The  increase  in  the 
magnitude  of  the  pressure  fluctuations  in  going  from  upstream 
towards  the  rotor  in  the  axial  direction  is  clear,  as  is  the 
rapid  decrease  at  the  rotor  outlet.  Outside  the  rotor  exit 
plane  the  signal  becomes  flat  and  that  there  sure  two  blades 
is  no  longer  obvious.  Signals  for  K9.5  and  K1Q  can  be  com¬ 
pared  with  pictures  taken  from  the  oscilloscope.  Fig.  32,  on 
which  the  time-base  can  be  adjusted  to  provide  a  record  of 
more  than  one  blade  pair.  Some  differences  in  signals  from 
different  passages  are  evident  which  are  probably  due  to 
small  geometric  differences  in  the  hardware.  In  particular, 
this  is  certainly  the  case  for  the  rotor  leading  edge  plane 
at  K9  as  can  be  seen  from  Fig.  53.  The  signals  from  four 
different  blade  pairs  are  plotted,  and  differences  in  the 
maximum  and  minimum  values  are  clearly  evident.  This  is  due 
to  slight  differences  in  the  tip  gap  and  in  the  chords  of 
the  blades.  It  is  noted  that  the  tips  of  the  blades  pass 
approximately  across  the  center  of  transducer  K9,  on  the 
average . 
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The  largest  amplitudes  occur  for  the  transducer  K9.5, 
Fig.  52.  Consequently  this  location  was  chosen  to  locate  a 
pneumatic  static  pressure  probe  of  special  shape  which  was  to 
be  used  in  conjunction  with  the  iterative  data  reduction  pro¬ 
cedure  described  in  [34].  The  geometric  shape  of  the  probe 
follows  exactly  the  suggestions  given  in  [34].  The  probe 
was  manufactured  so  that  it  could  be  substituted  at  any  wall 
transducer  location. 

A  comparison  of  the  measurements  obtained  with  the 
standard  pneumatic  static  pressure  tap,  the  diameter  of 
which  was  1  mm,  with  those  obtained  with  the  special  design 
showed  that  there  was  no  detectable  difference.  However,  it 
is  suggested  that  for  higher  rotor  speeds  with  the  attendent 
increase  in  fluctuation  amplitudes  the  comparison  should  be 
repeated. 

Additionally,  signals  from  one  of  the  transducers 
were  acquired  with  different  numbers  of  repetitions  and  a 
least  squares  curve  fit  method  of  up  to  the  order  of  10  was 
used  to  approximate  each  set  of  results.  This  was  done  in 
an  attempt  to  overcome  the  problem  of  long  acquisition  times 
when  the  initial  software  [6]  was  in  use.  In  principle,  of 
course,  the  same  procedure  could  be  applied  to  A  and  B 
probe  signals.  In  Fig.  54  results  are  shown  for  5  different 
numbers  of  repetitions  with  a  curve  fit  of  ninth  order.  The 
original  signal,  which  is  plotted  for  Rep  ■  1  only,  appears 
as  noise  about  the  fitted  curve.  It  is  quite  remarkable 


that  results  for  different  numbers  of  repetitions  are  almost 
the  same.  The  procedure  which  is  a  form  of  data  smoothing 
or  averaging  could  be  improved  either  by  dividing  the  results 
for  the  256  counts  into  more  than  2  sections  or  by  using  a 
more  sophisticated  mathematical  method  (which  might  result  in 
capacity  problems  for  the  21-MX  computer) .  Although  the 
software  improvements  [10]  have  reduced  test  time  remarkably, 
it  may  be  worthwhile  to  keep  this  procedure  in  mind  for  fu¬ 
ture  measurements  if  a  speed-up  of  the  data  acquisition  is 
required . 


VI.  PROPOSAL  FOR  A  FIRST  APPLICATION  IN  THE  RGG 
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The  test  facility  for  Rotating  Annular  Cascades  (RGG)  of 
the  DFVLR  was  first  used  for  investigation  on  cylindrical 
(conical  flow  angle  a  =  0°)  shaped  test-wheels  [1  and  2], 

In  this  configuration  the  test-section,  basically,  consisted 
of  two  cylindrical  rings  forming  an  annulus.  Following  this 
program,  since  1980  experiments  have  been  underway  dealing 
with  conical  cascade  flows.  Such  flow  conditions  occur,  for 
example,  in  the  last  stages  of  steam  turbines,  particularly 
in  the  tip  regions.  The  goal  was  to  obtain  information  on 
the  flow  behavior  for  those  rotor  areas  where  it  is  obvious 
that  flow  particles  can  not  travel  on  streamlines  which  are 
nearly  parallel  to  the  machine  axis  (which  is  the  commonly 
used  simplification  in  flow  models) . 

Three  new  test  sections  have  been  built  for  three  dif¬ 
ferent  conical  flow  conditions:  a  ■  30°  (hub  and  tip) , 
a  »  45°  (hub  and  tip)  and  a  *  40°  (hub)  and  a  »  45°  (tip)  . 

A  sketch  of  the  test  section  for  the  second  case  (a  *  45°) 
is  given  in  Fig.  55.  It  is  noted  that  measurement  planes 
and  probe  settings  are  chosen  in  the  manner  indicated  so  as 
to  be  able  to  make  full  use  of  available  hardware.  Details 
of  the  complex  program  which  is  carried  out  in  cooperation 
with  Forschungsvereinigung  Verbrennungskraftmaschinen,  e.  V. 
(FW) ,  Frankfurt,  West  Germany  are  given  (for  example)  in 
(47]. 
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Basically,  the  measurement  and  evaluation  method  used  for 
conical  flows  is  the  same  one  that  was  developed  for  cylin¬ 
drical  flows.  Of  course,  somewhat  different  assumptions  are 
involved  in  the  data  reduction  procedures  for  each  case.  In 
the  downstream  flow  field  only  those  quantities  which  are 
easiest  to  measure  need  be  evaluated;  namely,  absolute  total 
pressure  and  absolute  total  temperature.  Assuming  that 
there  exists  a  core-flow  area  in  the  axial  or  conical  direc¬ 
tion  of  the  annulus  where  only  minor  changes  in  flow  proper¬ 
ties  occur,  two  dimensional  cascade  data  can  be  obtained. 
Results  which  were  necessary  to  prove  this  assumption  for 
the  cylindrical  flow  conditions  are  given  in  [1]  and  [2]. 

In  [48]  it  is  demonstrated  that  for  conical  flow  conditions 
(a  »  30°)  such  a  core-flow  also  exists.  Nevertheless,  it 
would  be  highly  desirable  to  have  information  of  the  local 
time-dependent  flow  quantities  also. 

Among  several  motivations  the  most  important  is,  of 
course,  that  it  would  be  of  more  than  only  academic  interest 
to  verify  that  the  core-flow  areas  are  fundamentally  two  di¬ 
mensional  in  character,  and  to  use  an  independent,  second 
measurement  method  to  do  it.  Although  in  [1]  and  [2]  the 
limits  of  this  core-flow  were  defined  with  respect  to  blade 
shapes,  cascade  geometries  and  downstream  Mach  numbers,  it 
is  suggested  that  for  different  cascades  (=  test  wheels)  a 
check  of  those  regions  should  be  done.  Of  course,  it  is  not 
necessary  to  measure  point  by  point  throughout  the  different 
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flow  fields  of  interest,  but  some  limited  investigations 
should  be  carried  out.  It  is  expected  that  these  areas  can 
be  evaluated  both  more  accurately  and  more  quickly  (with  re¬ 
spect  to  test  time)  with  high  response  pressure  probes*. 
Therefore,  it  is  suggested  that  they  should  be  used  in  paral¬ 
lel  with  the  existing  probes  and  techniques  in  the  future. 

In  addition,  conditions  can  occur  where  the  region  of 
nearly  two-dimensional  flow  is  quite  small.  This  has  hap¬ 
pened  so  far  for  extreme  flow  conditions  (high  supersonic 
relative  downstream  Mich  numbers,  in  particular!  in  the 
cylindrical  case  as  well  as  in  the  conical  one.  Results 
carried  out  up  to  now  indicate  that  the  quality  of  two- 
dimensional  flow  decreases  with  increasing  conical  flow 
angle,  o.  The  investigations  for  o  »  45°  are  now  being 
prepared,  and  it  can  be  expected  that  the  flow  behavior 
will  continue  this  trend.  This  provides  a  further  reason 
why  preparations  should  be  made  to  make  available  a  second 
measurement  method  which  has  the  potential  of  getting  more 
detailed  information  on  the  flow. 

Finally,  it  may  also  be  possible  to  carry  out  measure¬ 
ments  of  general  interest  concerning  the  problem  of  how  to 
measure  time-averaged  pneumatic  pressures.  The  systematic 
investigations  which  are  available  so  far  are  limited  to 
fluctuation  frequencies  of  about  4  kHz.  This  is  due  to  the 
relatively  small  devices  which  have  been  used  to  produce  a 
periodically  interrupted  flow  field.  The  magnitude  of  the 
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tangential  velocities  is  limiting  and  the  need  to  have  the 
wave  length  equal  to  several  times  the  probe  dimension  has 
placed  a  practical  limit  on  the  available  frequency.  It 
would  be  worthwhile  to  extend  this  region  to  ranges  of 
about  10  kHz,  which  are  typical  characteristic  blade  pas¬ 
sing  frequencies  in  many  machines.  For  this  purpose  wheels 
in  the  test  section  of  the  RGG  could  serve  as  nonsteady 
flow  generators,  only.  By  using  different  kinds  of  probes, 
with  a  calibrated  high  response  probe  as  a  reference  device, 
the  comparison  of  results  for  the  same  axial  and  radial  lo¬ 
cations  (of  course,  in  different  circumferential  positions) 
can  be  expected  to  give  quantitative  results  of  interest. 
Another  advantage  of  the  test-section  is  that  in  the  down¬ 
stream  flow  field,  measurements  can  be  made  in  measurement 
planes  which  cover  a  wide  range  in  the  axial  direction. 

This  can  not  be  done  in  the  TX  without  removing  the  stator 
and  could  not  be  done  in  many  other  similar  research  fa¬ 
cilities.  Doing  similar  measurements  at  different  axial 
distances  could  lead  to  additional  information  concerning 
the  problem  of  time  averaging  fluctuating  pressures. 

In  these  investigations  it  is  suggested  that  an  oil- 
filled  measurement  device  such  as  is  described  in  [33]  and 
[34]  be  included.  The  reported  probe  is  well  understood 
and  proven  by  its  originators  who  have  agreed  to  demonstrate 
its  application  in  the  RGG. 


In  the  previous  section ,  estimates  were  made  with  the 
aid  of  examples  of  how  an  incorrectly  evaluated  time-averaged 
total  pressure  can  influence  results  obtained  using  the  DPDS 
technique.  Open  questions  are  evident ,  and  so  far  no  other 
experience  exists  concerning  the  behavior  of  the  A  and  B 
probes  in  areas  of  high  subsonic  and  supersonic  nonsteady 
f lows . 

The  following  remarks  concern  the  application  of  the 
DPDS  technique  in  the  RGGs  The  first  goal  should  be  to 
provide  additional  qualitative  information  of  the  complex 
conical  flow  field,  if  possible  for  the  current  tests  and 
if  this  can  not  be  done  in  the  time  available,  for  future 
experiments . 

The  blade  passing  frequencies  in  the  RGG  are  of  the  same 
order  of  magnitude  as  those  in  the  TX.  The  maximum  wheel 
speed  is  only  half  (250  revs/sec) ,  but  the  number  of  blades 
can  be  2.5  times  larger.  This  can  occur  for  a  tip-section 
cascade  with  a  pitch- to-chord  ratio  of  1.,  low  stagger  angle 
(about  20°)  and  low  incidence  angle  (of  the  same  range) . 

This  would  result  in  a  frequency  of  f  »  12.5  kHz,  which  can 
be  regarded  as  the  highest  frequency  to  be  investigated. 
Assuming  f  *  10  kHz  and  50  blades  the  maximum  speed  will  be 
200  revs/sec  (12000  RPM)  which  may  be  a  more  realistic  value. 

For  hub  section  cascades,  where  lower  pitch- to-chord 
ratios  are  used  (about  .5}  the  number  of  blades  is  about  100 
and,  again,  assuming  £  *  10  kHz  will  lead  to  speeds  of  100 


revs/sec  (6000  RPM) ,  which  is  about  double  what  has  been 
found  to  be  necessary  from  past  experience. 

First  of  all  it  is  necessary  to  provide  a  trigger  system 
for  the  data  acquisition.  In  principal,  a  variety  of  solu¬ 
tions  sure  possible.  A  sketch  in  Fig .  5 6  demonstrates  the 
one  which  is  suggested  for  the  RGG.  The  test-wheels  them¬ 
selves  act  as  timing  disks  and  two  photo-electrical  sensors 
(NANO-SCAN,  Skan-A-Matic  Corp. ,  Elbridge,  NY)  are  used  to 
provide  pulses  at  one-per-blade  and  one-per-revolution  fre¬ 
quencies.  The  installation  of  the  sensors  within  the  inner 
casing  (position  (X)  ,  Fig.  56)  is  expected  to  be  easy,  mainly 
because  the  sensors  have  a  flexible  "snout"  which  in  addi¬ 
tion  can  provide  an  adjustment  of  the  distance  between  the 
sensors  and  the  wheel.  For  initial  tests,  which  will  be 
carried  out  without  a  device  similar  or  equal  to  the  "Pacer," 
only  the  "one-per-revolution"  signal  is  needed. 

On  the  back  side  of  the  test-wheels  (see  Fig.  56)  radial 
straight  lines  are  required  to  be  coincident  with  the  cen¬ 
ters  of  the  blades.  This  must  be  done  accurately  to  provide 
a  stable  input  to  the  pacer  and  permit  data  acquisition  be¬ 
hind  blade  pairs  of  interest  independent  of  minor  wheel 
speed  variations.  Both  one-per-blade  and  one-per-rev.  time- 
mark-radii  should  be  chosen  to  be  as  large  as  possible,  to 
provide  the  most  distinct  electrical  pulses. 

The  complete  arrangement  can  be  checked  out  using  a  disk 
which  represents  the  wheel  and  speeding  it  up  to  test 
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conditions  (up  to  12,0Q0/min)  outside  the  test  section. 

Such  tests  would  provide  answers  concerning  the  optimum  dis¬ 
tance  between  sensor  and  wheel  (a  first  estimate  can  be  taken 
from  data  given  in  Skan-A-Matic’s  handbook)  and  will  also 
determine  what  width  of  the  marks  in  circumferential  direc¬ 
tion  is  necessary.  If  the  width  is  small  in  comparison 
to  the  distance  between  two  adjacent  blades  a  "wave  shaper" 
section  in  a  pacer  device  (see  Section  III)  may  not  be  neces¬ 
sary.  Consequently,  the  marks  should  first  of  all  be  made 
by  simply  laying  it  out  with  a  metal  scribe.  It  is  possible 
that  reasonable  results  could  be  obtained  using  a  small 
groove  filled  with  some  kind  of  reflecting  material. 

It  may  be  possible  to  use  the  blades  themselves  to  pro¬ 
vide  the  time-marks.  In  this  case  the  second  sensor  would 
be  installed  in  position  (b) ,  Fig.  56  and  from  position  (a) 
only  the  "one  per  revolution"  pulse  would  be  needed.  This 
procedure  should  be  used  if  it  is  more  complicated  than  ex¬ 
pected  to  provide  well  defined  time-marks  for  the  case 
as  described.  Complications  for  case  (b)  may  result  from 
the  abrasive  material  used  in  the  outer  casing  groove  to 
guarantee  optimum  tip  clearance.  Perhaps  particles  of  this 
material  would  disable  the  sensor.  In  any  event,  a  "wave- 
shaper”  would  be  necessary  because  the  geometry  of  the  blades 
would  produce  some  kind  of  a  oell-shaped  signal. 

The  necessary  investment  is  approximately  4000. -DM  which 
includes  two  units  and  a  third  one  for  replacement  purposes. 


In  Fig.  57  a  sketch  of  a  total  pressure  probe  which  can 
be  used  as  an  equivalent  to  the  A  probe  is  shown.  Although 
it  is  desirable  (or  may  be  necessary)  that  the  flow  should 
act  directly  on  the  head  of  the  sensor  and  even  small  vol¬ 
umes  should  be  avoided  because  of  resonance  problems,  geo¬ 
metric  restrictions  of  the  hardware  do  not  allow  such  a 
design  at  the  moment.  In  particular,  the  diameter  of  the 
holes  in  the  casing  where  the  probes  can  be  inserted  into 
the  flow  field  is  6mm  and,  of  course,  the  maximum  value  of 
the  axial  extension  of  a  probe  must  be  slightly  less  than 
this  value.  A  rearrangement  is  not  possible  at  the  present 
time.  Therefore,  the  proposed  modified  A  probe  is  a  compro¬ 
mise  between  what  is  possible  and  what  is  needed,  however  it 
is  expected  that  it  can  serve  adequately  for  familiarization 
with  the  new  technique. 

As  shown  in  Fig.  57,  the  probe  will  rotate  about  an  axis 
through  the  probe  tip.  Because  of  the  design,  manufacturing 
problems  might  occur,  particularly  since  the  transducers 
must  be  handled  extremely  carefully.  In  the  knee-area  (Fig. 
57)  a  hole  (for  the  transducer)  and  a  slot  (for  the  refer¬ 
ence  tube)  are  necessary  to  adjust  the  sensor  in  the  indi¬ 
cated  position.  An  adhesive  secures  the  transducer  in  that 
location  and  is  used  for  closing  the  hole  and  the  slot.  It 
may  be  worthwhile  to  use  dummy  transducers  initially  to  find 
out  whether  the  installation  will  work  in  the  suggested  way. 
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In  Fig.  58  a  second  probe  for  the  RGG  is  sketched  which 
can  serve  for  measurements  of  the  static  pressure  (see  for 
example  Fig.  32) .  First  of  all  this  device  should  be  used  to 
obtain  wall-static  pressures  in  the  different  measurement 
planes  where  probes  can  be  inserted.  In  addition,  radial 
probe  traverses  in  the  annulus  are  possible  using  the  probe 
because  the  dimensions  are  chosen  to  fit  the  traverse  mechan¬ 
isms.  This  may  be  most  worthwhile  in  the  upstream  regions 
where  the  pitch  angle,  <J>,  is  expected  to  be  small  but,  never¬ 
theless,  it  may  also  be  of  interest  to  try  it  downstream. 

In  comparison  to  a  transducer  which  is  flush-mounted  in 
the  usual  way  at  the  casing  wall  the  probe  in  Fig.  58  takes 
more  time  to  build  and  the  additional  expense  is  justified 
only  because  of  its  more  varied  use.  It  is,  however,  easier 
to  build  than  the  total  pressure  probe  shown  in  Fig.  57  and 
does  have  the  advantage  that  the  flow  acts  directly  on  the 
sensor. 

In  addition,  a  pneumatic  equivalent  to  the  static  pres¬ 
sure  probe  of  Fig.  58  should  be  manufactured.  This  will 
allow  an  on-line  calibration  as  discussed  in  Section  IV, 
and  comparisons  can  be  made  with  a  static  pressure  probe 
which  is  already  available.  In  the  annulus  downstream,  a 
five-hole  cone  probe  (giving  time  averaged  total  pressure, 
yaw  angle  and  pitch  angle)  would  provide  a  better  under¬ 
standing  of  the  time  averaged  flow  field. 
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The  RGG  can  also  be  used  for  the  calibration  of  probes. 

In  [1],  [2],  respectively  the  calibration  of  a  temperature 
probe  was  carried  out.  In  this  case,  the  test-wheel  was 
replaced  by  calibration  rings  which  produced  a  contour  over 
the  inner  cylindrical  wall.  The  calibration  rings  were  cal¬ 
culated  as  two-dimensional  half-nozzles.  Each  ring  was 
instrumented  along  the  contoured  surface  with  tappings  to 
obtain  surface  pressure  measurements  to  compare  with  the 
design  distributions.  Calibrations  in  the  subsonic  regime 
were  carried  out  using  a  single  ring  producing  a  simple 
contraction.  Different  Mach  numbers  were  produced  by  ad¬ 
justing  the  flow  rate  generated  by  the  radial  supply  com¬ 
pressor.  For  two  different  supersonic  Mach  numbers,  1.2 
and  1.4,  two  different  rings  were  constructed  to  designs 
which  were  produced  by  the  method  of  characteristics. 

In  the  case  of  conical  test-sections,  calibration  rings 
are  not  necessary.  The  test-wheel  is  to  be  replaced  by  a 
disk  which  serves  to  complete  the  inner  contour.  With  dif¬ 
ferent  compressor  speeds  and  different  bypass  ratios  the 
Mach  number  can  be  varied  because  the  annulus  is  acting  as 
a  nozzle.  In  the  measurement  planes  @  and  (l)  ,  Fig.  55, 
even  supersonic  speeds  can  be  set  as  can  be  seen  from  static 
pressure  measurements  along  the  inner  and  outer  contours 
[47]. 

The  Reynolds  number  can  be  controlled  independently  of 
the  Macy  number.  This  is  effected  by  controlling  independently 
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the  total  pressures  and  total  temperatures  in  the  closed 
cycle.  The  total  pressure  is  regulated  by  either  sucking 
with  a  rotary  compressor  of  sliding  vane  type  or  by  supplying 
air  through  the  bypass-duct.  The  temperature  is  controlled 
by  the  cooler  following  the  radial  supply  compressor  [1  and 
21. 

With  the  devices  and  test  procedures  described  in  the 
previous  paragraphs  it  should  be  possible  to  obtain  informa¬ 
tion  on  the  wake-flow  behavior  and  boundary  layer  effects  in 
the  inner  and  outer  casing  wall  areas.  Information  could 
be  obtained  initially  simply  by  recording  signals  photo¬ 
graphically  from  the  screen  of  an  oscilloscope,  as  demon¬ 
strated  in  Fig.  32. 

So  far,  no  provisions  are  discussed  for  the  evaluation 
of  pitch  angles,  <fr.  In  [48]  calculations  based  on  a  3- 
dimensional  time  marching  method  are  reported  which  indicate 
a  remarkable  streamline  shift  even  when  the  blade  height  is 
relatively  small.  In  the  example  given  the  flow  was  de¬ 
flected  from  *  45°  at  the  inlet  plane  to  a2  -  30°  down¬ 
stream.  These  numbers  may  be  used  to  estimate  the  a-ranges 
which  may  possibly  occur  in  the  experimental  conditions. 
However,  it  is  recommended  that  the  techniques  described 
abrve  be  applied  first,  and  a  device  to  enable  pitch  angle 
measurements  should  be  designed  later  when  a  background  of 
experience  in  the  RGG  has  been  established. 


VII.  CONCLUSION 


The  Dual  Probe  Digital  Sampling  (DPDS)  Technique  for 
measuring  three  dimensional  flow  fields  in  turbomachines  has 
been  described  and  its  development  at  the  Turbopropulsion 
Laboratory  in  recent  years  has  been  reviewed.  Analysis  of 
the  first  results  in  the  transonic  compressor  at  the  Naval 
Postgraduate  School  was  reported. 

Two  different  data  reduction  procedures — an  approximate 
analytical  and  a  numerical  one— were  examined.  When  tested 
on  data  obtained  in  free  jet  experiments,  the  numerical 
method,  which  used  calibration  data  arrays  rather  than 
analytical  approximations  of  calibration  data,  was  found  to 
give  the  better  results.  This  was  explained  as  being  the 
result  of  the  inability  of  fourth  order  polynomial  functions 
used  in  the  analytical  procedure  to  adequately  represent  the 
data  over  the  required  range  of  yaw  angles.  Nevertheless, 
the  application  of  the  numerical  method  of  reduction  to  data 
from  the  compressor  did  not  allow  the  velocity  vector  dis¬ 
tribution  to  be  determined  with  acceptable  certainty. 

In  the  analysis  of  the  initial  results,  which  is  reported 
here,  the  behavior  of  high  response  probes  as  a  function  of 
yaw  angle  in  a  steady  flow  free  jet  and  in  the  unsteady  flow 
behind  the  rotor  were  shown  to  be  different.  The  differences 
were  greatest  for  data  acquired  within  the  clearly  unsteady 
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region  of  the  rotor  blade  wakes.  Differences  were  least  for 
data  acquired  within  the  nearly  periodic  region  outside  the 
blade  wakes.  Consequently  it  was  concluded  that  measurements 
obtained  in  the  compressor  can  not  be  reduced  at  the  present 
time  by  a  straightforward  application  of  the  reported  methods 
to  the  acquired  data  using  the  calibration  carried  out  in 
steady  flow;  to  do  so  will  lead  to  incorrect  results  with  an 
unknown  uncertainty.  This  result  is  thought  to  have  impor¬ 
tant  implications  concerning  the  interpretation  of  unsteady 
data  from  more  complex  multiple-sensor  probe  geometries. 

Speculative  reasons  for  the  differences  in  behavior  were 
offered;  however,  time  was  not  available  to  the  author  to 
investigate  further.  The  full  understanding  must  be  left  ?s 
the  goal  of  the  future  program.  It  is  suggested  that  close 
collaboration  between  and  coordination  of  the  programs  at 
the  Naval  Postgraduate  School  and  at  the  DFVLR-AVA  Gottingen 
Fluid  Dynamics  Branch,  where  similar  measurements  are  to  be 
made  in  rotating  cascades,  will  lead  most  efficiently  to  the 
required  understanding  and  to  further  development  of  the  DPDS 
and  other  high  response  probing  techniques. 
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Figure  5 :  Configuration  of  probes  and  locations  of  case-wall  transducers 
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Figure  12:  Frequency,  generated  with  wave-generator;  f  -  10  kH 
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Results  of  analytical  data  reduction  applied 
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Figure  52:  Case-wall  results;  Rep  *  20 
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